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Abstract 

The  classical  engineering  theory  of  bending  due  to  Bernoulli  and 
Euler  serves  as  a  cornerstone  for  structural  analysis  and  design. 
Limitations  of  this  theory,  however  become  apparent  in  flexural  wave 
propagation  studies;  it  predicts  infinite  phase  velocity  as  the  wavelength 
becomes  shorter.  This  theoretical  deficiency  is  corrected  by  Timoshenko 
theory  which  accounts  for  transverse  shear  deformation.  A  thorough  study 
of  several  exact  elasticity  solutions  reveals  that  there  are  two  additional 
effects  that  are  of  the  same  order  as  transverse  shear  in  bending  behavior. 
These  are  due  to  (a)  transverse  normal  strain  and  (b)  an  additional 
cfohrr  ibut  ion  to”  axial  s Ere AS.  \  A  new  engineering  bending  theory  which 
-eeecmn£e:H?«r---the*e  is  presented.  Predictions  of  static  bending  response 
using  this  theory  agree  exactly  with  elasticity  solutions  for  several 
uniformly  distributed  loading  cases.  The  contributions  due  to  various 
physical  effects  are  found  to  be  more  pronounced  for  orthotropic  materials 
with  low  shear  and  transverse  extensional  moduli.  Such  properties  are 
typical  of  advanced  composite  materials  used  in  the  aerospace  industry. 
The  theory  is  extended  to  study  the  dynamic  behavior  of  beams  and  static 
buckling  of  columns. 

Validation  for  the  theory  is  provided  by  analysis  of  classic 
benchmark  problems  -  -  -  a  simply  supported  beam  under  sinusoidally 
distributed  loading  and  flexural  wave  propagation  in  rectangular  slabs. 
Numerous  applications  are  presented.  Effects  of  property  degradation  due 
to  hygrothermal  conditioning  on  the  behavior  of  several  elementary 
unidirectional  composite  structures  are  studied.  The  hygrothermal 
condition  used  to  simulate  long  term  aircraft  service  does  not  pose  a 
serious  problem.  The  loss  of  performance  is  approximately  ten  percent. 
However,  complete  moisture  saturation  produces  significant  effects. 
Reductions  in  the  performance  up  to  35  percent  can  be  expected,  thus 
indicating  that  this  condition  should  be  avoided.  - 
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SUMMARY 

The  classical  engineering  theory  of  bending  due  to  Bernoulli  and 
Euler  dates  back  to  1705  and  precedes  the  theory  of  elasticity  by  over  100 
years.  It  has  long  been  recognized  as  a  convenient  approximation  for 
slender  beams  and  serves  as  a  cornerstone  for  structural  analysis  and 
design. 

Limitations  of  engineering  bending  theory  become  apparent  in 
studying  the  propagation  of  elastic  flexural  waves  of  short  wavelength, 
i  The  Bernoulli-Euler  theory  predicts  infinite  phase  velocity  for  harmonic 

waves  as  the  wavelength  becomes  shorter.  This  result  is  of  course, 
physically  absurd.  This  theoretical  deficiency  is  corrected  by  the  theory 
j  proposed  by  Timoshenko.  In  Timoshenko  theory,  the  influence  of  transverse 

shear  deformations  are  accounted  for,  which  results  in  a  finite  limit  for 
phase  velocity. 

|  A  thorough  study  of  several  exact  elasticity  solutions  reveals  that 

there  are  two  additional  effects  that  are  of  the  same  order  as  transverse 
shear  in  bending  behavior.  These  are  due  to  transverse  normal  strain  and 
t  an  additional  term  in  the  axial  stress.  A  new  engineering  theory  of  planar 

bending  «Aiich  accounts  for  these  is  presented  in  this  work.  Predictions  of 
static  beam  bending  response  using  the  new  equations  agree  exactly  with 
<  elasticity  solutions  for  several  uniformly  distributed  loading  cases.  The 

theory  is  validated  by  means  of  a  thorough  consistency  analysis  and  by 
comparing  with  an  exact  solution  to  a  nonuniform  loading  case. 


The  theory  ie  extended  to  dynamics  and  validated  through  a 
consistency  analysis.  It  is  applied  to  study  flexural  wave  propagation  in 
slabs  and  vibration  behavior  of  beams.  The  theory  is  validated 
quantitatively  by  comparing  with  the  classic  benchmark  problem  -  -  - 
flexural  wave  propagation  in  slabs.  The  results  indicated  superior  range 
of  applicability  compared  to  Timoshenko  or  Bernoulli-Euler  theories.  The 
theory  is  further  established  by  extracting  Stephen  and  Levinson's  theory 
specialised  to  thin  rectangular  beams  from  the  new  dynamic  equations. 

An  elementary  theory  to  static  buckling  analysis  and  preliminary 
estimates  of  buckling  loads  are  provided  for  simply  supported  columns. 
Results  are  in  general  agreement  with  Timoshenko  theoretical  predictions. 

Applications  of  practical  interest  are  provided  through  the  study 
of  hygrothermal  effects  on  the  flexural  behavior  of  composite  beams  of 
unidirectional  layup.  Hygrothermal ly  degraded  mechanical  properties  are 
used  in  cosqmting  the  response  under  static  and  dynamic  situations. 

A  sumary  of  the  conclusions  based  on  the  results  and  suggestions 
for  future  work  are  provided. 


CHAPTER  1 


INTRODUCTION 

Use  of  fiber  reinforced  resin  matrix  composite  materials  in 
aerospace  vehicles  is  increasing.  This  is  primarily  due  to  their  superior 
mechanical  properties  and  the  ease  with  which  they  can  be  tailored  to  a 
specific  application.  The  properties  of  composites  depend  on  the 
individual  properties  of  the  constituents  and  the  manner  in  which  the 
fibers  are  utilized.  The  most  structurally  efficient  type  of  laminated 
composite  is  cosq>osed  of  layers  of  unidirectional  continuous  fibers.  In 
this  case,  mechanical  properties  depend  also  upon  the  fiber  orientation, 
which  may  be  chosen  arbitrarily.  This  permits  tailoring  to  specific  design 
requirements. 

The  directional  nature  of  the  composite  material's  mechanical 
properties  poses  unique  challenges  for  the  analyst.  Consider,  for  example, 
a  single  layer  or  lamina  made  of  a  composite  material.  The  extensional 
modulus  along  the  direction  of  fibers  is  usually  very  large  relative  to  the 
extensional  moduli  in  the  lateral  directions  and  the  shear  moduli.  This  is 
a  marked  departure  from  conventional  isotropic  materials.  The  result  is 
that  the  relative  importance  of  physical  effects  is  influenced  by  the 
directional  nature  of  properties  and  their  relative  magnitude.  Transverse 
shear  deformations,  for  exaaple,  are  much  more  pronounced  for  composite 
structures. 

Transverse  shear  deformation  effects  in  connection  with  beam  and 
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plate  bending  have  been  studied  extensively.  However,  there  still  is  no 
unique  way  of  accounting  for  them.  An  engineering  theory  which  includes 
them  in  a  simple,  rational  way  is  desirable.  This  is  a  primary  objective 
of  the  present  work. 

Currently,  considerable  research  activity  in  the  area  of  composite 
materials  is  directed  towards  the  study  of  hygrothermal  effects  and  three- 
dimensional  effects  such  as  delamination.  Resin  matrix  materials  absorb 
moisture,  particularly  in  elevated  temperature  environments.  As  a  result, 
the  matrix  softens  and  matrix  controlled  properties  show  significant 
degradation.  This  is  due  to  a  lowering  of  glass  transition  temperature  of 
the  resin  matrix  material.  The  resulting  degradation  of  stif fness-related 
and  strength-related  properties  is  a  serious  problem  for  designers. 

Other  problems  of  considerable  concern  are  attributed  to  three- 
dimensional  effects.  Two  such  problems  are  matrix  micro-cracking  and 
delaaination.  Analytical  solutions  are  accomplished  by  using  three- 
dimensional  numerical  techniques.  These  solutions  are  very  expensive  to 
construct  and  are  often  inaccurate  in  transition  regions.  Interlaminar 
shear  stresses  and  transverse  normal  stresses  are  thought  to  be  the  primary 
causes  of  the  aforementioned  failures.  Several  theories  have  been  proposed 
recently  to  determine  these  stresses  more  accurately.  The  resulting 
equations  are  cumbersome  and  the  results  are  not  fully  satisfactory.  For  a 
preliminary  design  analysis,  an  engineering  theory  that  is  simple  yet 
reliable  would  be  a  positive  contribution. 

An  historical  discussion  of  bending  theory  is  presented  to 
establish  the  basis  for  new  developments  and  to  permit  the  present  work  to 
be  placed  in  proper  perspective.  Second,  an  analysis  of  bending  behavior 
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is  described  which  utilizes  an  exact  solution  from  the  theory  of  elasticity 
for  isotropic  materials.  A  unique  feature  is  the  use  of  tracer  constants 
in  order  to  track  the  contributions  due  to  various  physical  effects 
throughout  the  course  of  the  analysis.  With  the  aid  of  insight  from  this 
analysis,  a  new  engineering  theory  is  proposed.  The  theory  is  applied  to 
elementary  static  applications  for  beam-type  structures  which  illustrate 
its  use  and  permit  comparisons  with  exact  elasticity  solutions  to  establish 
its  validity.  It  is  extended  to  study  the  dynamic  behavior  of  beams  and 
static  buckling  of  columns. 

The  theory  is  applied  to  study  the  effects  of  the  property 
degradation  due  to  hygrothermal  conditioning  on  composite  structural 
behavior  under  static  and  dynamic  loading  situations. 

The  scope  of  this  work  is  restricted  to  planar  bending  situations. 
In  its  present  form,  the  theory  applies  to  beams  with  thin  rectangular 
cross  sections  which  respond  to  planar  bending  in  plane  stress  or  to 
infinitely  wide  plates  which  respond  in  plane  strain  (cylindrical 
bending).  Both  isotropic  and  orthotropic  materials  are  considered.  Beams 
of  orthotropic  material  are  the  simplest  type  of  structures  where  composite 
material  behavior  can  be  studied. 


CHAPTER  II 


HISTORICAL  SKETCH 

I 


Introductory  Remarks 

A  brief  history  of  the  development  of  bending  theory  is  given  below. 
Although  the  emphasis  is  on  engineering-type  bending  theories,  some 
studies  involving  three-dimensional  exact  elasticity  equations  and  higher 
order  beam  and  plate  bending  theories  are  included. 
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Classical  Theories  of  Bending 

The  detailed  historical  development  of  the  mathematical  theory  of 

1  2 

elasticity  is  given  in  the  books  by  Love  ,  Todhunter  and  Pearson  ,  and 
3 

Sokolnikoff  .  The  classical  theory  of  planar  bending  of  beams  is  due  to 
Bernoulli  and  Euler.  James  Bernoulli  derived  the  relationship  between 
bending  moment  and  curvature  in  1705.  Euler  assumed  this  relation  in  his 
analysis  of  the  elastica  and  vibration  of  thin  rods  .  However,  the  full 
engineering  bending  theory  in  its  present  form  is  due  to  Coulonfo  .  Coulomb 
clarified  the  equilibrium  equations  and  introduced  the  notion  of  neutral 
axis.  The  theory  is  based  on  the  hypothesis  that  plane  sections  normal  to 
the  neutral  axis  remain  plane  and  unextended  after  bending.  It  provides  a 
convenient  approximation  for  slender  beams  and  serves  as  a  cornerstone  for 
structural  analysis  and  design. 

The  classical  engineering  theory  of  plate  bending  had  its  origin  in 
the  pioneering  work  of  Sophie  Germain.  She  was  awarded  a  prize  in  1815  for 
her  attempt  to  provide  a  theoretical  basis  for  the  modal  figures  obtained 
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in  Chladni's  vibration  experiments^.  Her  work  was  finally  published  in 
4 

1821  .  It  contained  an  error  in  the  expression  for  strain  energy  of 

2 

bending,  which  was  corrected  by  Lagrange  .  The  governing  differential 

equation  for  flexural  vibration  of  plates  was  independently  established  by 

5  6  7  8 

Navier  ,  Poisson  ,  and  Cauchy  .  However,  it  was  Kirchhoft  who  resolved 

the  famous  controversy  concerning  the  nature  and  number  of  proper  boundary 
conditions.  Love*  provided  an  extension  for  the  bending  of  shells. 

The  Kirchhoff-Love  theory  of  plate  and  shell  bending  is  based  on  the 
hypothesis  that  normals  to  the  neutral  surface  remain  normal  and 
unstretched  after  bending.  This  permits  only  two  boundary  conditions  per 
edge.  Three  boundary  conditions  per  edge,  however,  provide  a  more 
realistic  behavioral  description  of  the  plate  and  shell  bending. 

Limitations  of  elementary  bending  theory  become  apparent  in  stu¬ 
dying  the  propagation  of  elastic  waves  of  short  wavelengths.  It  was 

9 

pointed  out  by  Lamb  that  Bernoulli-Euler  theory  is  inadequate  for  impact 
type  loads.  It  leads  to  the  physically  absurd  conclusion  that  disturbances 
are  propagated  instantaneously  throughout  the  beam.  This  is  because  it 
predicts  infinite  phase  velocity  for  harmonic  waves  as  the  wavelength 
becomes  shorter.  According  to  the  exact  solution  of  Rayleigh*®,  this 
should  approach  a  finite  limit.  Rayleigh  attempted  to  improve  the 
classical  beam  bending  theory  by  accounting  for  rotatory  inertia  effects 
and  obtained  a  finite  limit. 


Elasticity  Solutions 


The  adequacy  of  a  specific  theory  can  best  be  decided  by  comparing 
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it  with  exact  solutions.  Solutions  of  the  full  three-dimensional 

elasticity  equations  are  rare.  Usually  simple,  closed  form  solutions  do 

not  exist.  Some  can  be  found,  however,  for  a  few  sufficiently  simple 

geometric  configurations  and  for  simple  loadings.  Normally,  some 

numerical  technique  is  required  to  solve  a  practical  problem. 

Reference  1  gives  some  of  the  earliest  solutions  for  beams  bent 

under  arbitrarily  continuous  loading.  Closed  form  solutions  for  beams  of 

arbitrary  cross-section  bent  by  terminal  couples  and  loads  are  provided  by 
2  11  12 

Saint  Venant  .  Pochhamner  and  Chree  have  studied  independently  wave 

propagation  in  an  infinitely  long  beam  of  solid  circular  cross-section.  A 

similar  study  for  an  infinitely  wide  rectangular  plate  was  done  by  Lord 

Rayleigh^.  Hia  solution  for  the  flexural  wave  velocity  of  propagation  in 

an  isotropic  rectangular  slab  is  a  bench  mark. 

13  14 

Von  Karman  and  Seewald  undertook  studies  of  the  flexure  of 
rectangular  beams.  Their  attempts  to  correlate  elasticity  solutions  with 
the  classical  beam  bending  theory  yield  corrections  to  the  moment- 
curvature  relationship.  It  appears  that  Pearson*  was  first  to  report  these 
corrections.  Later,  Grashof,  Michell  and  Filon*  also  provided  similar 
corrections  for  beams  subjected  to  distributed  loading  independently. 
These  are  later  attributed  to  shear  deformation  by  several  analysts. 
Goodier*^,  however,  in  his  three-dimensional  order  of  magnitude  analysis 
of  beam  bending,  showed  that  the  correction  term  is  not  necessarily  due 
only  to  transverse  shear  stress  effects.  His  analysis  indicates  that  the 
transverse  normal  stresses  and  an  additional  term  in  the  expression  for 
axial  stress  may  also  contribute  to  beam  bending  response.  No  means  of 


accounting  for  these  effects  is  offered,  however.  It  is  rather  perplexing 
that  later  researchers  have  not  followed  up  on  Goodier's  work  or  made  an 
effort  account  for  the  aforementioned  effects.  A  recent  exception  is 
Reference  36. 

Using  a  different  approach,  Donnell1**’17  obtained  solutions  like 

18 

those  of  References  13  and  14.  Hashin  proposed  a  simple  direct  method  to 

obtain  compatible  stress  field  in  beams  subjected  to  polynomial  loading. 

Later,  he  extended  it  to  obtain  exact  stresses  in  plane  orthotropic 
19 

beams  .  Hashin1 s  approach  differs  from  earlier  work  in  that  there  is  no 

guess  work  involved  in  obtaining  solutions.  In  the  earlier  work,  the 

solutions  are  guessed  or  found  by  combining  known  solutions  so  as  to 

satisfy  boundary  conditions.  Though  a  direct  method  was  first  proposed  by 
20 

Neou  ,  the  choice  of  the  degree  of  polynomial  remained  arbitrary  in  his 
study . 

21 

Recently  Cheng  has  provided  a  plate  theory  based  upon  the  three- 
dimensional  equations  of  elasticity.  He  did  not  consider  any  transverse 
loading  in  the  development. 

Shear  Deformation  Theories 

Attempts  to  remove  the  theoretical  shortcomings  in  classical 
bending  theories  gave  rise  to  theories  incorporating  certain  refinements. 
During  vibration,  beam  cross  sections  experience  rotatory  motions  as  well 
as  translations.  Also,  the  transverse  deflection  of  a  beam  has 
contributions  due  to  the  transverse  shearing  forces  as  well  as  the  bending 

moment.  Correction  for  the  influence  of  rotatory  inertia  was  provided  by 

10  22  22 
Rayleigh  as  mentioned  earlier.  Grashof  (1878)  and  Rankine  (1895) 

included  the  effects  of  transverse  shear  deformation  in  analysing  some 

static  beam  bending  problems. 


A  more  refined  theory  that  accounts  for  rotatory  inertia  and 

transverse  shear  deformation  was  proposed  by  Timoshenko  in  his  famous 
23 

paper  in  1921.  Timoshenko's  theory  is  widely  recognized  and  used 
wherever  improvement  on  classical  theory  is  sought.  It  is  of  some  his¬ 
torical  interest  that  both  the  rotatory  inertia  correction  and  transverse 
shear  correction  were  given  first  by  the  French  analyst  M.  Bresse  in  1859 
in  his  Cours  de  Mecanique  Appliquee.  This  work  has  been  overlooked  in  the 
later  development  of  the  subject. 

In  Timoshenko's  theory,  there  is  a  shear  correction  factor,  k,  to 

account  for  nonuniform  shear  stresses  across  the  cross-section  of  the  beam. 

Originally  it  was  taken  as  the  ratio  of  average  to  maximum  shear  stress  on 

the  cross-section.  Thus,  for  a  rectangular  cross-section,  this  procedure 

24 

yields  a  value  of  2/3  for  k.  In  a  subsequent  paper  ,  Timoshenko  proposed 
a  new  value,  8/9,  for  better  correlation  with  the  experimental  results  of 
Filon^. 

Most  of  the  refined  theories  that  followed  are  based  upon  Timo- 
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shenko-type  beam  equations.  Uflyand4  and  Mindlin  have  developed  plate 

theories  including  the  effects  of  rotatory  inertia  and  transverse  shear 

deformation  specifically  for  dynamic  applications.  Uflyand's  equations 

are  essentially  an  extension  of  Timoskenko's  beam  equations. 

26 

Mindlin' 8  plate  theory  also  contains  the  shear  correction  factor, 
k.  A  unique  way  of  obtaining  k  is  not  provided  in  the  theory.  It  is  chosen 
by  an  ad-hoc  criterion.  He  suggested  two  such  criteria  based  on  a  matching 
principle.  One  is  to  choose  k  such  that  the  limiting  phase  velocity  for 
very  short  flexural  waves  is  made  identical  with  the  velocity  of  Rayleigh's 
surface  waves.  The  other  is  to  select  k  so  that  it  gives  exact  circular 


frequency  of  the  first  antisymmetric  mode  of  thickness-shear  vibration. 

The  result  is  that  k  depends  on  the  cross  sectional  shape  and  the  mode  of 

motion.  Values  of  k  for  various  cross  sectional  shapes  are  provided  in 

Reference  27.  References  28  and  29  treat  the  analysis  of  thickness  shear 

vibration  in  quartz  crystals  where  Kindlin' s  shear  deformation  theory  was 
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used.  Traill-Nash,  and  Collar  and  Goodman  and  Sutherland  have  provided 
analogous  theories  for  beam  vibration.  In  the  former,  experimental 
verification  was  also  given. 

Following  the  above  mentioned  work,  there  were  several  attempts  to 

improvise  Mindlin-type  equations.  The  differences  primarily  relate  to  the 

selection  of  shear  correction  factor,  k,  according  to  various  ad-hoc 

32 

criteria.  A  recent  survey  by  Kaneko  gives  an  excellent  description  of 
33 

them.  Cowper  provided  a  new  formula  for  k  which  depends  upon  Poisson's 

ratio.  His  equations  appeared  to  give  satisfactory  results  for  static 

applications  and  long-wavelength,  low  frequency  deformation  of  beams.  In 

the  analysis,  the  effects  of  transverse  normal  stress  were  neglected. 
.  .  34 

Leibowitz  and  Kennard  have  used  an  alternative  approach  to  obtain  k. 
Exact  bending  moment-curvature  relation  of  a  beam  bent  under  its  own  weight 
was  used  to  redefine  k. 

A  shear  deformation  theory  with  two  arbitrary  constants  is 
presented  in  Reference  35.  The  underlying  idea  was  further  developed  in 
Reference  36  by  combining  the  contributions  of  References  33,  34  and  35. 
Two  constants,  k^  and  kj*  are  introduced  and  chosen  as  follows,  kj  is  the 
usual  shear  correction  factor  defined  by  a  transverse  shear  stress-strain 
relationship.  k.  is  associated  with  bending  stiffness  and  comes  from 


accounting  for  transverse  normal  and  lateral  stresses  in  the  moment- 

curvature  relation.  It  is  assumed  that  the  stresses  during  beam  vibration 

can  be  approximated  by  those  of  a  beam  bent  under  uniform  gravity  loading. 

36 

The  above  work  of  Stephen  and  Levinson  is  especially  noteworthy. 
It  is  the  first  work  to  make  use  of  the  observations  and  insights  of 
Goodier^.  Good  agreement  with  exact  solutions  for  several  flexural  wave 
propagation  problems  has  been  demonstrated  for  overall  response 

properties.  The  theory  requires  an  exact  elasticity  solution  for  St. 
Venant's  bending  problem  for  the  beam  cross  section  under  consideration  for 
its  application,  however.  It  also  possesses  a  shortcoming.  The 

displacement  variables  are  averaged  quantities  over  the  cross  section. 

Consequently,  no  claims  are  made  regarding  the  pointwise  distribution  of 
stresses  or  displacements. 

37 

For  static  applications,  Reissner  has  derived  plate  ^quitiont 

that  account  for  shear  deformation  and  transverse  normal  stresses  using  an 

entirely  different  approach  in  1945.  Consequently,  his  shear  deformation 

theory  is  marked  by  the  presence  of  transverse  normal  stress  effects. 

He  used  assumed  stresses  and  an  energy  principle  to  obtain  the  governing 

38 

differential  equations.  He  clarified  and  subsequently  improved  his 

39 

earlier  work  for  bending  of  plates  without  transverse  normal  stress 

40 

Several  analogous  theories  for  plates  are  surveyed  in  Fane’s  book  .  The 
notable  among  them  are  due  to  Pane,  Hencky  and  Kroon. 

Refined  Bending  Theories 

The  classical  bending  theory  and  the  shear  deformation  theories  can 
be  considered  to  be  related.  The  former  can  be  obtained  by  taking  the 


shear  modulus  in  terms  associated  with  the  transverse  shear  deformation  to 
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be  very  large  in  the  latter.  In  both  theories,  the  displacement  variables 
contain  only  linear  terms  in  the  thickness  coordinate,  z.  More  refined 
theories  can  be  formulated  by  taking  higher  order  terms  in  z  in  addition  to 
the  linear  terms.  These  theories  are  sometimes  called  higher  order 
theories. 

The  usual  procedure  for  developing  this  type  of  theory  has  been  a 
displacement  formulation.  The  displacements  are  assumed  as  a  power  series 
in  the  thickness  coordinate,  z.  The  governing  differential  equations  and 
the  appropriate  boundary  conditions  are  then  obtained  using  energy 
principles.  References  41  and  42  formed  the  basis  for  subsequent  works. 
Whitney  and  Sun^  and  Nelson  and  Lorch^  have  proposed  refined  theories  for 
laminated  plates  and  shell  structures.  A  class  of  contact  problems  in 
beams  is  solved  in  Reference  45.  An  excellent  survey  of  various  higher 
order  theories  and  a  comparative  study  of  relative  differences  is  given  in 
Reference  46.  Lo,  Christensen  and  Wu  ^  have  proposed  recently  a  theory 
for  isotropic  and  laminated  composite  plates  specifically  for  dealing  with 
problems  rfiich  involve  rapidly  fluctuating  loads  with  a  characteristic 
length  of  the  order  of  thickness. 

Displacement  formulations  of  the  above  types  begin  with 
kinematically  admissible  displacements,  but  the  stress  equilibrium 
equations  are  violated.  In  Reference  48,  a  systematic  approach  to  obtain 
solutions  of  the  three-dimensional  elasticity  equations  is  given  for  beams 
subjected  to  arbitrary  loading.  It  is  generally  true,  however,  that 
assumed  displacement  approaches  result  in  stresses  which  are  not  in 
equilibrium  and  which  provide  poor  design  stress  estimates.  To  obtain 
stresses,  alternative  means  may  have  to  be  sought  as  pointed  out  in 


Reference  49. 
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Reissner's  *  ’  plate  theory  and  the  theory  proposed  in  the 

Reference  36  should  also  be  classified  as  refined  theories  of  bending. 
This  is  because  effects  besides  those  due  to  transverse  shear  deformation 
are  accounted  for. 


Although  engineering  bending  theory  has  a  long  and  successful 
history,  there  are  issues  that  remain  surrounded  in  uncertainty.  Beyond 
the  classical  theories,  there  are  no  unique  or  clearly  superior  refined  or 
shear  deformation  theories.  A  thorough  exploration  of  the  effects 
enumerated  by  Goodier^  in  addition  to  those  of  transverse  shear  has  not 
been  undertaken.  Are  there  relationships  among  the  various  theories  that 
have  been  proposed  and  are  some  "better"  than  others?  Problems  and 
opportunities  associated  with  composite  structures  require  that  these 


issues  be  resolved. 


CHAPTER  III 


PRELIMINARY  ANALYSIS 


Introductory  Remarks 

As  a  first,  important  step  in  an  analysis  of  bending  behavior,  a 
plane  stress  elasticity  solution  for  a  simply  supported  beam  under  uniform¬ 
ly  distributed  loading  is  studied.  It  is  possible  to  identify  the  in¬ 
dividual  contributions  due  to  various  factors  affecting  beam  response.  An 
assessment  of  their  relative  importance,  therefore,  can  be  made. 

Problem  Definition  and  Solution 

The  two-dimensional  elasticity  solution  for  a  simply  supported  beam 

under  uniformly  distributed  loading  ia  given  in  the  text  by  Timoshenko  and 
22 

Goodier  .  It  is  valid  for  very  thin  rectangular  beams  in  the  plane  stress 
form.  For  infinitely  wide  plates,  the  same  solution  remains  valid  if  a 
transformation  of  elastic  constants  for  plane  strain  is  employed. 

The  beam  and  coordinate  system  are  shown  in  Figure  l.  The  length  of 
the  beam  is  21,  and  the  depth  is  2c.  The  width  of  the  beam  is  taken  as  unity 
for  convenience.  The  beam  is  bent  by  a  uniformly  distributed  load  of 
intensity  q  applied  to  its  upper  surface.  The  midspan  of  the  beam 
centroidal  axis  is  chosen  as  the  origin  for  the  coordinate  axes  x  and  z.  z 
»  4c  and  z  ■  -c  correspond  to  the  bottom  and  top  surfaces  of  the  beam.  The 
notation  and  convention  are  shown. 

For  the  stresses,  the  usual  convention  and  notation  are  followed. 
Accordingly,  o  ia  the  axial  stress,  0  is  the  transverse  normal  stress 
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ando^  ia  the  transverse  shear  stress.  They  are  given  by  the  following 
expressions : 


y«»  -  ft  <*2  -  *  §i  <^r  -  ^)  in 

°zz  "  “  fa  (f*  "  °2z  +  T“  >  (2) 

~  §1  z2^x  (3) 

3 

I  is  the  second  moment  of  the  cross  sectional  area  and  is  2c  /3  for  the 

rectangular  section  under  consideration.  These  stresses  satisfy  all  the 

governing  differential  equations  and  the  stress  boundary  conditions  on  the 

upper  and  lower  surfaces.  On  the  ends  x  *  H,  the  stress  boundary 

conditions  are  satisfied  in  an  overall  Saint  Venant  sense. 

As  an  aid  in  this  analysis ,  three  tracer  constants,  a  ,a  and  a  , 

an  s 

are  introduced.  They  are  defined  and  used  so  as  to  facilitate  keeping 
track  of  three  distinct  contributions  to  the  response.  The  first  term  of 
Equation  (l)  corresponds  to  the  bending  stress  giveu  by  classical 
Bernoulli-Euler  theory.  The  underlined  term  is  a  stress  contribution  which 
will  be  called  the  "nonclasaical  axial  stress".  It  produces  no  resultant 
force  or  moment  and  is,  therefore,  a  self  equilibrating  stress,  a  is  the 

d 

tracer  constant  associated  with  this  contribution.  If  a  *  1,  this 
contribution  is  fully  accounted  for.  If  *  0  in  the  following,  however, 
it  is  ignored  and  the  Bernoulli  -  Euler  axial  stress  distribution  is 
recovered.  For  example,  the  axial  stress  is  written  using  this  convention 
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a  n  and  ag  are  defined  analogously  and  are  associated  with  contributions  due 

to  o  and  a  ,  respectively, 
zz  XZ  r 

The  displacement  components  u  and  w  are  shown  in  Figure  l.  Expres¬ 
sions  for  them  can  be  obtained  by  using  Hooke’s  law  and  the  strain  dis¬ 
placement  relations.  For  plane  stress,  Hooke's  law  for  an  isotropic 
material  is 
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E  is  Young' 8  modulus,  \>  is  Poisson's  ratio  and  G  is  the  shear  modulus.  z  xx 
and  e  are  the  extensional  strains  in  the  longitudinal  and  the  transverse 
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directions  and  y  is  the  shear  strain.  The  strain-displacement  relations 
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With  the  aid  of  Equations  (1) 
integration  of  Equations  (7). 
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*1 ,  -  i  are  imposed : 


(6)  ,  u  and  w  can  be  obtained  by  direct 
The  following  boundary  conditions  at  the 
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w(i,0)  *  w(-l,0)  ■  0 


(8) 


They  represent  support  conditions  applied  at  the  beam  axis.  Also,  from  the 
symmetry  requirement 


u(0,z)  »  0 


(9) 


These  conditions  are  sufficient  to  prevent  rigid  body  motion. 

The  expressions  for  u  and  w  are 

3  2 

u(x,z)  *  £(|2X  -  y-)z  +  aa(^  2)x 

3  3 

+  anV(L3'  c 2*  +T-)x]  (10) 

w(x,z)  *  w(x,0)  -  J|y[%  -  +  an 

7  7  —2  4  2  2  7 

♦ * (r  -  x2)  ;  +  wa(|-  -  (in 


In  Equation  (11),  w(x,0)  is  the  vertical  deflection  of  the  beam  centroidal 
axis  due  to  bending.  It  is  given  by  the  Equation 


w(x,0)  “8  -  -  YJ  +  [d+,,)as  “  +  ~y~)]c2x2J  (12 

as 


inhere 
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8  is  the  deflection  at  the  midspan  of  the  beam. 


An  Analysis  of  Beam  Response 

|  The  major  differences  between  the  classical  Bernoulli-Euler  theory 

and  the  elasticity  solution  can  be  clearly  identified  in  Equation  (13). 

The  underlined  term  represents  the  correction  to  the  former  due  to  the 

|  presence  of  contributions  identified  by  the  tracer  constants  a^,  anda^. 

Note  that  the  contributions  due  to  all  three  effects  -  transverse  shear, 

nonclassical  axial  stress  and  transverse  normal  strain  -  are  of  the  same 

23 

'  order  of  magnitude.  A  static  version  of  Timoshenko's  theory  includes 

only  the  terms  associated  with  a  . 

1  a 

The  corrections  shown  in  Equation  (13)  were  known  to  earlier 
13  14 

authors  ’  .  However,  they  did  not  differentiate  among  the  various 

contributions.  This  differentiation  provides  the  key  ingredient  for  the 

establishment  of  a  rational  engineering  theory.  Goodier^  suspected  that 

|  the  other  influences  beside  transverse  shear  were  important,  but  offered  no 

means  of  estimating  them  quantitatively  and  no  concrete  examples  of  their 

contribution  to  beam  response.  The  approach  adopted  here  makes  the  matter 

<  transparent  and  settles  the  issue  for  this  example. 

If  V  is  taken  to  be  0.3  and  CX  =  a  -a  =>  1  in  Equation  (13),  then 

n  8  a  n 
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(14) 


A  corresponding  result  from  Timoshenko's  original  shear  deformation 
23 

theory  can  be  obtained  by  setting  an  “  aa  *  and  a8  =  1  in  Equation 
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(13).  For  V  =  0.3,  the  result  is 

^Timoshenko  "  [l+  3‘12  J7 ]  °5) 

4 

A  popular  alternative  is  to  use  Timoshenko  theory  with  c*b  =  -j.  This 
approximately  corresponds  to  result  obtained  if  Reissner's  approach  is 
adopted.  This  leads  to 

^  Timoshenko  "  WeY  ['  *  fl]  <154> 

Conclusions 

On  the  basis  of  the  foregoing  analysis,  the  following  conclusions 
are  reached: 

1.  A  Timoshenko-type  transverse  shear  theory  does  not  contain  the 
necessary  physical  ingredients  to  treat  problems  with  distributed 
loadings. 

2.  Transverse  shear,  nonclassical  axial  stress  and  transverse 
normal  strain  make  contributions  to  the  response  that  are  of  the  same  order 
of  magnitude.  A  theory  that  is  purported  to  be  more  accurate  or  complete 
than  classical  theory  must,  therefore,  correctly  account  for  all  of  these 


influences . 


CHAPTER  IV 


FOUNDATIONS  OF  A  NEW  THEORY 

Objectives 

The  primary  abjective  of  this  work  is  the  development  of  a 
foundation  for  an  engineering  bending  theory  which  is  consistent, 
reliable,  and  simple  to  use.  The  theory  should  provide  more  reliable 
information  than  existing  ones.  Furthermore,  it  must  account  for  the  three 

effects  that  were  clearly  identified  previously  -  transverse  shear 

strain,  nonclassical  axial  stress  and  transverse  normal  strain. 

An  engineering  theory  is  one  in  Which  assumptions  or  approximations 
are  introduced  in  order  to  simplify  the  governing  equations  or  facilitate 
their  solution.  Hopefully  only  a  little  accuracy  is  sacrificed  for  a 
considerable  reduction  in  computational  labor.  The  intent  is  to  encompass 
the  heart  of  the  problem  under  consideration.  Consistency  and  rationality 
are  desirable,  but  mathematical  rigor  is  meaningless  in  this  context.  An 
engineering  theory  is  judged  solely  on  the  basis  of  the  results  obtained 
from  its  use. 

The  standard  of  comparison  for  results  that  is  used  herein  is 
rigorous  solution  to  the  equations  of  elasticity  theory  for  the  problem  in 
question. 

A  second  objective  is  to  obtain  stress  estimates  that  are 
improvements  over  those  provided  by  classical  bending  theory.  This  must  be 
accomplished  if  the  influence  of  nonclassical  axial  stress  is  to  be 
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properly  accounted  for. 


Statically  Equivalent  Stresses 

Equilibrium  of  a  beam  element  is  governed  by  overall  equations 
containing  resultant  axial  force,  shear  force  and  bending  moment .  The  sign 
convention  and  notation  for  these  appear  in  Figure  2.  The  equilibrium 
equations  are 
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The  force  and  moment  resultants  are  defined  in  terms  of  stresses  as 
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In  the  above,  a  rectangular  cross  section  of  unit  width  is  assumed  as 
before.  In  addition,  the  beam  is  assumed  to  be  of  uniform  depth. 

*  According  to  classical  theory,  the  stresses  are 
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A  is  Che  cross  sectional  area,  which  is  2c  for  Che  rectangular  cross 
section  under  consideration.  These  stresses  are  statically  equivalent  to 
the  applied  loads  and  satisfy  the  stress  equations  of  equilibrium: 


a  +o  *  o 
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In  addition,  Equations  (19)-(21)  are  satisfied,  as  are  appropriate 
stress  conditions  at  z  ■  c  and  z  ■  -c. 

The  above  stresses,  although  not  exact,  serve  as  a  first  approxi¬ 
mation.  This  stress  field  is  statically  equivalent  to  applied  loads, 
however,  it  does  not  satisfy  compatibility  requirements.  It  will  be  used 
subsequently  to  develop  approximations  for  the  displacement  components. 


Kinematics 

Classical  Bernoulli-Euler  theory  is  based  upon  a  kinematic  as¬ 
sumption  that  is  equivalent  to  ignoring,  and  hence  setting  to  zero, 
transverse  normal  strain  and  transverse  shear  strain.  Timoshenko-type 
shear  de format' on  theories  account  for  transverse  shear  strain  but  still  do 
not  permit  transverse  normal  strain.  On  the  basis  of  the  previous  analysis 
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of  the  simply  supported  beam  example,  it  appears  necessary  to  completely 
abandon  the  Bernoulli-Euler  kinematic  assumption.  In  order  to  obtain  some 
simplification  from  the  complete  elasticity  equations,  an  assumption  that 
facilitates  the  analysis  is  required,  however. 

The  central  assumption  that  replaces  the  Bernoulli-Euler  hypothesis 
in  the  present  development  is  that  the  statically  equivalent  stresses  in 
Equations  (22)-(24)  can  be  used  to  estimate  the  transverse  normal  strain 
and  transverse  shear  strain.  Note  that  this  is  an  assumption  regarding 
stresses.  It  is  not  a  kinematic  assumption.  This  is  in  sharp  contrast  to 
classical  and  Timoshenko-type  shear  deformation  theories. 

The  development  will  be  carried  out  for  orthotropic  materials  with 
principal  material  directions  corresponding  to  axes  of  the  beam.  The 
appropriate  form  of  Hooke's  Law  for  plane  stress  (beams  of  thin  rectangular 
c ross  section)  is 
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C  ,  e  ,  and  e  are  the  axial  strain,  transverse  normal  strain  and 
XX  ’  zz’  xz 

transverse  shear  strain,  respectively.  and  E^  are  elastic  moduli 

associated  with  the  x  and  z  directions.  V  ^  is  Poisson's  ratio  and  is 


the  transverse  shear  modulus. 


Although  it  is  impossible  to  obtain  a  unique  displacement  field  from 
the  incompatible  stress  field,  a  selective  use  of  strain  displacement 
relations  and  Hooke's  Law  permits  an  approximate  form  for  the  displacement 
field  to  be  determined.  The  error  involved  in  this  process  will  be 
estimated  in  Chapter  VII. 

On  the  basis  of  the  above,  Equations  (7),  (22),  (24)  and  (28)  permit 
the  transverse  normal  strain  to  be  approximated  as 
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Integration  of  this  equation  results  in  the  following  expression  for  the 
lateral  displacement  component,  w: 


w  *»  W(x) 
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(31) 


W(x)  is  the  lateral  deflection  of  the  beam  axis  (*  »  0),  which  is  an  unknown 
function  to  be  determined. 

The  axial  component  of  displacement  u  can  be  estimated  as  follows. 
Equations  (7),  (16),  (18),  (23),  (29),  and  (31)  permit  u  to  be  expressed 

t  * 


in  terms  of  the  shear  stress  and  w 
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This  expression  is  integrated  to  yield 


u 


U(x)-2W 


»X 


v  3 
.  13  n  z 

+  £71  Q  T 


2G 


13 


,  2 

t  (c  2  -  r 


2  3 

c  zJ  1  3  2, 

— r-~  +  -^  c  z  ) 


(32) 


U(x)  is  the  axial  deflection  of  the  beam  axis,  which  is  an  unknown  function 
to  be  determined. 

The  static  displacement  field  is  completely  described  by  Equations 
(31)  and  (32).  U  and  W,  the  axis  displacement  components,  emerge  as 
natural  kinematic  variables.  If  and  '*,QO  in  (31)  and  (32),  a 

tranverse  shear  theory  is  obtained  which  includes  the  effects  of  cross 
section  warping.  If,  in  addition,  Gj^>00>  then  the  classical  Bernoulli- 
Euler  kinematic  assumption  is  recovered. 

Considerable  simplification  is  achieved  if  the  underlined  terms  in 
Equations  (31)  and  (32)  are  neglected.  These  terms  are  associated  with 
higher  derivatives  of  the  shear  force  Q  than  the  remaining  terms.  This 
simplification  is  adopted  here.  Its  full  implication  will  be  discussed  in 
Chapter  VII.  The  accuracy  of  this  approximation  is  related  to  how  rapidly 
the  applied  load  q  varies  with  x. 
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Refined  Axial  Stress  Distribution 

The  axial  stress  <7  is  the  largest  and  most  important  stress 
component.  An  accurate  knowledge  of  it  is  often  all  that  is  needed  in  a 
practical  application.  A  refined  estimate  which  improves  Equation  (22)  is 
central,  therefore,  to  the  improvements  that  are  sought. 

Equations  (7),  (24),  (27)  and  (32)  can  be  utilized  to  prc-'.uce  a 
refined  axial  stress  expression. 
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In  the  above,  contributions  due  to  the  underlined  terms  in  Equation  (32) 
are  not  included.  Notice  that  the  stress  distribution  throughout  the 
thickness  is  not  linear  as  in  the  classical  approximation  (22). 

Relationships  for  the  axial  force  and  bending  moment  are  obtained  by 
using  Equations  (19),  (21)  and  (33).  The  results  are 
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The  parameter  kx  is  (Ej^/2Gj^  -V^)»  un^ty  ^or  an  isotropic  ma¬ 

terial.  Equations  (34)  and  (35)  permit  (33)  to  be  rewritten  as 


a 


XX 


,3  2 

(y  c  z  - 


( 33A) 


The  underlined  term  is  the  nonclassical  axial  stress  contribution,  wnich  is 
the  desired  refinement. 


The  governing  equations  for  the  new  theory  can  be  summarized  now. 
They  encompass  four  categories.  Overall  beam-type  equations  consist  of  the 
equilibrium  equations  ( 16)— ( 18)  and  the  constitutive  equations  (34)  and 
(35).  In  addition,  two  sets  of  equations  provide  the  distributions  of 
stresses  and  displacements  throughout  the  structure.  The  first  set  for 
stresses  consists  of  Equations  (33A),  (23)  and  (24).  The  second  for 

displacements  is  composed  of  Equations  (31)  and  (32)  with  the  underlined 
terms  omitted. 

The  above  collection  of  equations  requires  the  specification  of 
boundary  conditions.  The  classical  boundary  condition  options  are  to 
specify  N  or  U,  Q  or  W,  and  M  or  <J>  at  the  ends  of  the  beam.  $  is  a  rotation- 
related  variable. 

Three  commonly  used  rotation-related  variables  are  considered 
below.  The  first  is  the  rotation  of  the  cross  section  at  the  beam, 
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Another  is  the  rotation-related  variable,  which  is  defined  by 
the  following  equation: 
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This  variable  naturally  arises  in  Reissner's  development  of  plate  bending 
38 

theory  based  upon  the  complementary  energy  principle. 

The  third  is  the  mean  rotation  of  the  cross  section,  <J> 
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In  a  Timoshenko-type  theory,  since  u  is  linear  in  z,  all  of  the 

above  definitions  are  equivalent.  Equation  (36)  is  the  actual  definition 
.  .  23 

used  in  the  original  paper.  These  variables  permit  different  models  for 
simulating  clamped  end  conditions  to  be  defined. 

Discussion 

The  development  of  the  equations  requires  no  ad  hoc  kinematic 
assumptions  or  use  of  a  variational  principle.  The  central  assumption  is 
that  the  transverse  normal  and  shear  strain  components  can  be  estimated 
from  classical  stresses.  A  selective  use  of  strain  displacement  relations 


is  utilized  to  establish  the  approximate  form  of  the  displacements.  The 
equations  have  the  following  properties: 

1.  Stress  and  displacement  distributions  throughout  the  structure 
are  found  in  terms  of  the  response  variables  associated  with 
the  axis; 

2.  Nonclassical  axial  stress  and  cross  section  warping  effects, 
transverse  shear  strain  and  transverse  normal  strain  are  all 
accounted  for  in  a  rational  manner; 

3.  The  equations  can  be  shown  to  yield  exact  results  for  the  case 
of  uniformly  distributed  lateral  loading; 

4.  For  nonuniform  loading,  some  of  the  equations  are  approximate 
—  the  stresses  are  not  exactly  in  equilibrium  and  the  stresses 
and  displacements  are  not  exactly  compatible;  and 

5.  The  equations  are  as  simple  to  apply  as  static  Timoshenko-type 
shear  deformation  theories. 

Items  1-3  and  5  are  strong  points  in  favor  of  the  new  equations.  Item  4 

imposes  some  limitations  on  the  validity  of  the  theory,  which  will  be 

thoroughly  discussed  in  Chapter  VII,  but  it  is  responsible  for  the 

simplicity  that  is  achieved.  The  level  of  stress  approximation  which 

14 

results  is  analogous  to  that  suggested  by  Seewald  for  isotropic 
materials. 

In  the  process  of  solving  a  particular  bending  problem,  the  only 
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apparent  difference  from  application  of  a  static  Timoshenko-type  shear 

deformation  theory  is  the  value  for  the  coefficient  of  the  Q  -  term  in 

»x 

Equation  (35).  As  the  applications  will  demonstrate,  this  seemingly  minor 
difference,  together  with  the  use  of  Equations  (31),  (32)  and  (33A), 

produces  significantly  improved  results. 
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CHAPTER  V 


STATIC  APPLICATIONS 

Introductory  Remarks 

In  order  to  illustrate  the  benefits  of  the  new  theory,  several 

elementary  applications  for  uniform  beams  subjected  to  uniformly 

distributed  loading  applied  to  the  upper  surface  analogous  to  the 

situation  shown  in  Figure  1,  will  be  presented.  One  special  case  of  a 

linearly  varying  load  is  studied  to  illustrate  a  particular  point. 

Comparisons  are  made  with  the  exact  elasticity  solution,  classical 

23 

Bernoulli-Euler  theory  and  the  original  Timoshenko  theory  in  each  case. 
The  two  dimensional  elasticity  solution,  for  the  present  purposes,  is 
considered  an  exact  solution,  although  the  plane  stress  approximation 
requires  the  width  to  depth  ratio  of  the  beam  to  be  small.  A  discussion  of 
this  issue  is  given  in  Reference  23,  page  274. 

Solutions  are  derived  for  orthotropic  beams,  and  corresponding 
results  for  isotropic  beams  are  obtained  by  specialization.  Poisson's 
ratio  is  taken  to  be  0.3  throughout.  For  orthotropic  beams,  E  ^  i  13 
taken  to  be  30;  this  is  a  typical  value  for  a  modern  graph i te/ epoxy 
composite  material. 

In  presenting  results,  appropriate  response  variables  are  non- 
dimensionalized  with  respect  to  the  corresponding  values  obtained  from 
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Bernoulli-Euler  theory.  This  practice  permits  easy  recognition  of 
departures  from  classical  theory  predictions. 

Response  can  be  separated  into  bending  and  stretching.  Stretching 
is  governed  by  Equations  (16)  and  (34),  bending  by  Equations  (17),  (18), 
and  (35).  The  bending  problem  involving  M,  Q  and  W  must  be  solved  first.  N 
and  U,  stretching  variables,  are  determined  secondarily.  For  the  present 
purposes,  only  the  bending  portion  of  the  response  is  discussed. 

Simply  Supported  Beam 

The  exact  solution  for  a  simply  supported  (SS)  isotropic  beam  was 
presented  earlier.  The  precise  boundary  conditions  that  have  been  imposed 
at  the  ends  are 


SS:  M  -  0,  W  -  0 


(39) 


Such  a  beam  is  shown  in  Figure  1.  It  is  a  statically  determinate  struc¬ 
ture,  so  the  moment  and  shear  distributions  are  known. 

The  response  of  the  beam  is  defined  if  the  axis  lateral  deflection 
W(x)  is  found.  For  this  type  of  end  restraint,  W  can  be  expressed  as 
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In  the  above,  the  constant  K  is 


(41) 


The  tracer  constants  introduced  earlier  are  utilized  to  identify  the  origin 
of  the  various  contributions  to  K.  6  is  the  maximum  or  midspan  deflection, 
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which  is 
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L  *  2  is  Che  Cotal  lengCh  of  Che  beam  and  H  *  2c  is  Che  depCh  of  Che  cross 
seccion.  The  scress  distribuCion  aC  midspan  is  given  by 


(40B) 


The  soluCion  by  Che  presenC  cheory  correspond  Co  ot  =  ot  *  *  1  in 

(41);  iC  is  exacC  for  chis  problem.  If  K  is  seC  to  zero  in  Equations  (40) 

and  (42).  then  the  Bernoulli-Euler  result  is  obtained.  If  ot  *=  a  =0  and 

an 

«  *  1  in  (41),  the  staCic  Timoshenko  theory  prediction  is  recovered.  Note 

3 

that  Timoshenko  theory  overestimates  the  midspan  deflection  in  this  case. 


Cantilever  Beam 

For  a  cantilever  beam,  it  is  convenient  to  take  the  origin  of 
coordinates,  x  ■  0,  at  the  free  end.  x  ■  L,  corresponds  to  the  clamped  end. 
Unlike  the  more  elementary  theories,  the  present  theory  does  not  suggest  a 
unique,  simple  model  for  a  clamped  or  fixed  end.  Three  rotation  variables 
were  introduced  earlier  in  Equations  ( 36 ) — ( 38 ) .  Three  definitions  of 
clamping,  therefore,  will  be  discussed. 

All  results  can  be  cast  in  a  conmon  format.  The  three  types  of 
clamping  are  denoted  Cl,  C2,  and  C3.  They  correspond  to  the  following 
definitions: 
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Cl: 

W  = 

0, 

♦j  -  0 

(43) 

C2: 

W  = 

0, 

o 

II 

(44) 

C3: 

W  = 

0, 

*3  =  0 

(45) 

The  cantilever  beam  is  statically  determinate  with  M  and  Q  vanishing  at  the 
free  end.  At  the  fixed  end,  one  of  the  above  definitions  of  clamping  must 
be  imposed.  The  rotation  variables  can  be  expressed  in  the  common  form 
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The  constants  are 
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(46) 


(47) 


(48) 

(49) 


Since  the  structure  is  statically  determinate,  the  axial  stress 
distribution  can  be  readily  obtained  by  appropriate  substitution  for  N,  M 
and  Q  in  Equation  (33A).  The  result  for  the  fixed  end  is 
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(33B) 


The  lateral  deflection  can  be  conveniently  expressed  as  follows. 
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The  constants  5^  (i  =  1,  2,  3)  are  the  beam  tip  deflections;  they  are  found 
from 


5  . 
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[l  *  <2Ki  -  B  <f)2]  ; 


2,  3.  (51) 


The  constant  K  is  defined  in  Equation  (41). 

The  present  solutions  correspond  to  setting  all  the  tracer 
constants  to  unity;  they  are  exact  for  the  end  conditions  imposed.  The 
designation  Cl,  C2  and  C3  has  been  chosen  to  correspond  to  the  order  of 
increasing  stiffness  of  the  end  restraint.  The  coefficients  of  (H/L)  for 
Cl,  C2  and  C3  are  18.08,  12.15  and  8.19,  respectively,  for  the  orthotropic 
material  chosen. 

If  and  K  are  set  to  sero  in  the  above  equations,  the  Bernoulli- 

Euler  results  are  obtained.  This  approximation,  of  course,  overestimates 

stiffness.  The  Timoshenko  theory  result  is  obtained  from  the  Cl  case  by 

setting  a  *  1  and  a  “a  =*  0  in  Equations  (41)  and  (47).  The  coefficient 
s  an 

2 

of  (H/L)  is  15  for  the  same  material  considered  earlier.  It  is 
interesting  to  note  that  Timoshenko  theory  underestimates  maximum 
deflection  in  this  case. 

Bernoulli-Euler  theory  tends  to  always  overestimate  stiffness. 
Timoshenko  theory,  however,  in  light  of  the  results  presented  here,  may 
either  provide  an  overestimate  or  underestimate  of  the  maximum  deflection, 
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depending  upon  the  problem  under  consideration.  It  is  therefore,  "unre¬ 
liable"  in  this  sense. 

A  related  problem  is  a  cantilever  beam  subjected  to  a  linearly 
varying  distributed  load  that  varies  from  zero  at  the  free  end  to  q  at  the 
fixed  end.  The  exact  solution  is  given  in  Reference  22.  For  an  isotropic 
beam  with  Cl  restraint  at  the  fixed  end,  the  tip  deflection  is 
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(52) 


If  the  present  theory  is  used,  the  underlined  term  is  not  obtained.  For 
L/H  >2,  this  term  is  negligible.  For  practical  purposes,  therefore,  the 
present  theory  results  are  indistinguishable  from  the  exact  ones. 


Clamped  Beam 

Unlike  the  previous  examples,  the  clamped  beam  is  statically 
indeterminate.  Three  solutions  were  found  corresponding  to  the  three 
definitions  of  clamping  given  in  Equations  (43)-(45).  They  may  be 
expressed  in  a  common  form.  It  is  convenient  to  place  the  origin  of 
coordinates  at  midspan  and  use  the  semi-length  l  .  The  bending  moment 
distribution  is 


M  -  f  (|2  -  x2)  -  Mo 


(53) 


Mq  is  the  end  fixing  moment ,  which  is  positive  if  it  tends  to  reduce  the  end 


rotation  due  to  the  uniform  load 


The  expression  for  the  lateral  axis  deflection  is 


(54  A) 


j 

The  redundant  end  fixing  moment  is  different  for  each  type  of  clamping.  It 
can  be  written  in  the  conmon  form 
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(55) 


The  axial  stress  distribution  is  also  different  for  each  type  of  clamping. 
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As  before,  L  is  the  total  beam  length  and  Kj-K^  are  define<i  i°  Equations 
(47)-(49).  The  occurrence  of  different  end  moment  values  is  due  to  the 
statically  indeterminate  nature  of  the  structure.  The  present  theory 
yields  the  exact  solutions  to  this  problem  for  each  form  of  clamping. 

An  end  moment  ratio  as  a  function  of  beam  slenderness  is  plotted  in 
Figures  3  and  4  for  isotropic  and  orthotropic  materials,  respectively .The 
subscript  "B-E"  refers  to  the  value  from  Bernoulli-Euler  theory. 
Bernoulli-Euler  and  Timoshenko  theories  give  identical  predictions.  The 
present  theory,  however,  which  is  exact,  predicts  fundamentally  different 
behavior  that  differs  for  each  type  of  clamping.  Departures  from  classical 


theory  are  much  greater  for  beams  made  of  the  typical  orthotropic  material. 

An  interesting  phenomenon  occurs  for  Cl  restraint  of  orthotropic 
beams.  The  end  fixing  moment  acutally  reverses  sign  for  relatively  deep 
beams .  This  intriguing  situation  is  explained  by  the  fact  that  the  end 
rotation  is  forced  to  zero  by  a  combination  of  shear  force  and  bending 
moment.  By  virtue  of  synaetry  of  loading  and  structure,  the  end  shear 
force  is  fixed  by  vertical  force  equilibrium  considerations  alone. 
Consequently,  only  the  end  moment  is  available  for  controlling  rotation. 
Since  the  transverse  shear  stiffness  to  extensional  stiffness  ratio  is 
quite  low  for  this  material,  a  reversal  of  axment  is  required  to  offset  the 
large  shear  strain  at  the  axis  for  shorter  Cl-supported  beams.  A 
countertrend  for  C3~supported  beams  reflects  the  increased  relative 
difficulty  of  achieving  this  rigid  type  of  fixity  as  shorter  beams  are 
considered. 


Propped  Cantilever  Beam 

Let  the  origin  of  coordinates  be  the  simply  supported  end  of  a 
propped  cantilever  beam  and  x  *  L  be  the  clamped  end.  Three  cases 
corresponding  to  the  three  types  of  clamping  have  been  considered.  The 
bending  moment  is 


M ■  v " q  t 


Qo  is  the  shear  force  (reaction)  at  the  propped  end  x  ■  0.  It  can  be 
expressed  in  conon  form  as  follows  for  each  type  of  clamping. 
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The  lateral  axis  deflection  and  the  axial  stress  distribution  are  given  by 
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End  moment  ratio  plots  appear  in  Figures  5  and  6  for  isotropic  and 
orthotropic  beams,  respectively.  for  this  indeterminate  structural 
system,  Bernoulli-Euler  and  Timoshenko  theories  predict  different 
behavioral  trends.  The  present  theory  predictions  are  again  exact  for  this 
problem,  as  will  be  the  case  always  for  uniformly  distributed  loadings. 

Cl-supported  orthotropic  beams  again  exhibit  a  reversal  of  sign  of 
the  end  fixing  moment  similar  to  the  clamped  case.  Timoshenko  theory, 
which  approximates  this  end  fixity  condition,  displays  a  similar  trend,  but 
does  not  predict  an  actual  reversal  for  values  of  slenderness  parameter 
shown. 

Concluding  Remarks 

Several  representative  static  applications  whi:h  illustrate  the  use 
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of  the  new  equations  have  been  studied.  The  additional  effects  are  seen  to 
be  more  pronounced  for  statically  indeterminate  and  orthotropic 
structures.  Furthermore,  the  three  elementary  clamping  model  solutions 
indicate  that  care  must  be  devoted  to  matching  mathematical  descriptions  of 
boundary  restraint  with  practical  end  restraint  achieved  in  tests  or 
structural  assemblies.  The  sensitivity  of  the  response  to  boundary 
restraint  modeling  is  substantial  for  orthotropic  structures. 
Consequently,  the  next  chapter  is  devoted  to  a  study  of  this  issue. 
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CHAPTER  VI 

AN  APPROACH  TO  ACHIEVE  IDEAL  CLAMPING 


Preliminary  Remarks 

The  static  response  of  a  clamped  beam  under  uniformly  distributed 
loading  is  extremely  sensitive  to  the  precise  definition  of  the  boundary 
conditions.  This  sensitivity  is  more  pronounced  if  the  beam  is 
orthotropic.  None  of  the  three  elementary  models  for  the  clamped  end 
satisfy  exactly  the  generally  accepted  definition  of  zero  displacement  at 
the  fixed  end.  This  fact  is  illustrated  in  Figure  7.  It  shows  clamped  end 
cross  section  warping  for  a  typical  orthotropic  beam  with  L/H  “  4.  In  this 
Chapter,  an  approach  to  eliminate  the  warping  at  the  ends,  thereby 
achieving  ideal  clamping,  is  described. 


Analysis 

The  analysis  is  based  on  the  principle  of  superposition.  The  first 
part  of  the  solution  is  taken  to  be  one  of  the  elementary  clamping  models. 
To  this  solution  a  second  solution  for  the  beam  bent  by  prescribed  end 
displacements  is  added.  The  boundary  displacements  are  chosen  such  that 
they  nullify  the  warping  due  to  the  elementary  clamping  model.  Any  of  the 
three  elementary  models  may  be  chosen  as  the  starting  point.  However,  for 
the  purpose  of  illustration,  the  Cl  clamping  model  solution  is  taken  here 
as  a  starting  point.  Accordii  to  this  solution,  the  clamped  end  axial 
displacement  components  are 
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The  second  part  is  a  solution  to  the  following  boundary  value  problem: 


u  +  u1  =0 

at 

x  =  +  d 

(60) 

w  ♦  Wj  =0 

at 

x  =  +  d 

u  and  w  are  prescribed  boundary  displacements. 


0  -a  =0  on  z  =  +  c  for  all  x  (61) 

xz  zz  — 


The  addition  of  the  two  solutions  leads  to  a  refined  clamping  model  which 
will  be  referred  as  C4. 

ThiL  boundary  value  problem  is  solved  with  the  aid  of  the  principle 
of  virtual  complementary  work.  This  principle  is  appropriate  for  problems 
with  prescribed  displacement  boundary  conditions.  The  statement  is 
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The  principle  requires  stresses  which  satisfy  the  equilibrium  equations 


and  stress  boundary  conditions.  The  following  stress  field  is  selected. 
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B  is  an  arbitrary  constant  and  f2(x)  is  an  arbitrary  function  to  be 

determined  by  application  of  the  principle.  The  distribution  of  0  has 

been  selected  based  on  the  physical  nature  of  the  problem. 0  and  a  are 

xz  zz 

then  obtained  by  using  the  equilibrium  equations  (25)  and  (26)  and  the 

boundary  conditions  given  in  Equations  (61).  The  use  of  assumed  solutions 
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employing  free  functions  is  due  to  Kantorovich 

Substitution  of  Equations  (59)  and  (63)  into  (62)  leads  to  the 
following  functional: 
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Simplification  and  rearrangement  of  the  terms  of  Equation  (62A)  results  in 
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The  above  functional  yields  the  Euler  equations 
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where  yis  the  ratio  E^j/E^. 


The  natural  boundary  conditions  are 
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Equations  (65)-(68)  have  precise  physical  signif icance.  The  Euler 
Equations  are  the  relations  to  be  satisfied  for  the  kinematic  compatibility 
of  the  strains.  Equations  (67)-(68)  are  displacement  type  boundary 
conditions.  The  first  represents  transverse  shear  strain.  The  second 
represents  axial  strain. 
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T 

The  solution  to  Equation  (65)  may  be  chosen  as  e  ,  which  yields 


the  following  characteristic  equation 
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Consequently 
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The  roots  are  real,  equal  or  complex  depending  on  whether 


k*/y  >  2.045,  *  2.045,  <2.045 

i 

In  terms  of  elastic  constants  the  above  is 
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/(En/E33)>  2.045,  *  2.045,  <  2.045 


For  an  orthotropic  material  with  properties  Eu/Gi3  =  30, 

2 

E11/E33  =  15  and  v^3  =  0.3,  k^/  y  ■  14.406.  Consequently,  the  roots  are 
real  and  the  most  general  form  of  solution  to  Equation  (65)  may  be  written 
as 
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the  above  reduces  to 
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nij  and  are  to  be  obtained  from 
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For  an  isotropic  material  kx/y=l. 


The  roots  are 


4  o 

(+  3.656  —  +  1.528  — ).  The  solution  form  is.  therefore  selected  as 
—  c  —  c 
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f2  ■  C  cosh  —j—  cos  — j-  *  D  sinh  -j—  sin  —j- 


(71) 


where  and  are  given  by 
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(70A) 


It  is  interesting  to  note  that  for  isotropic  materials,  the  solution 
is  independent  of  the  material  constants  E,  G  and  V.  Also,  there  exists  a 
possibility  of  complex  and  equal  roots  for  orthotropic  materials  as 
indicated  by  Equation  (65B).  However,  this  situation  is  not  usually 
encountered  in  practical  situations.  The  practical  ranges  for  the 
parameters  kx  and  y  are  20-50  and  10-25,  respectively. 

Results  and  Discussion 

The  axial  stress  distribution  in  C4  restraint  beam  is 
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M  ^  is  the  end  moment  due  to  Cl  clamping.  From  Equation  (55)  it  may  be 


written  as 
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Equations  (71),  (72)  and  (66)  can  be  combined  to  cast  0  into  the 

xxC4 

following  convenient  form: 
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J  is  the  maximum  axial  stress  at  the  center  according  to  Bernoul ) i-Eul er 
theory  and  o  ^  i®  the  axial  stress  distribution  in  C2  restraint  beam. 
They  are  given  by 
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The  underlined  term  in  Equation  (73)  represents  a  correction  to  be  added  to 
C2  stress  distribution.  This  will  be  referred  as  local  disturbance 
parameter  in  the  subsequent  text. 
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It  is  possible  to  show,  for  large  values  of  L/H,  that  f^  can  be 
represented  approximately  by 

-»,C  _  -«JT 

f2  -  C  L/H  e  ♦  DL/H  e  (76) 

C  and  D  are  constants  and  £  is  the  dimensionless  distance  measured  from  the 
end. 

I  -  1  -  |  ( 76A) 

Equations  (73)  and  (76)  allow  the  following  to  be  written 

'  THp  <!>  {*V  '*»'  *  r} 
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LDP  refers  to  the  local  disturbance  parameter. 

An  important  conclusion  can  be  reached  by  observation  of  the 
Equation  (77).  The  local  disturbance  decays  exponentially  from  the  ends; 
for  relatively  slender  beams,  therefore,  the  solution  due  to  C4  approaches 
C2  results  in  the  interior  sone.  This  decay  phenomenon  is  well  understood 
and  is  usually  termed  as  an  end  effect  or  boundary  layer  effect^.  The 
maximum  value  of  the  local  disturbance  parameter  is  shown  graphically  in 
Figure  8  for  several  length-to-depth  ratios.  Results  are  obtained  from  the 
complete  expression  (69).  Another  point  of  interest  is  the  maximum  axial 
stress  at  the  edges  and  the  influence  of  the  orthotropicity  on  it.  Figure 
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9  provides  this  comparison.  In  the  calculations,  the  ratio  E^/E^-}  is 

taken  to  be  equal  to  E^/2G^  and  is  varied  between  10  to  50. 

The  isotropic  value  is  also  shown  for  comparison. 

The  presence  of  a  boundary  zone  near  ends  limits  validity  of  the 

elementary  theory  to  the  interior  zone.  The  lower  of  the  two  exponents 

and  m2  primarily  governs  the  extent  of  the  boundary  zone.  An  approximate 

estimate  of  a  decay  length,  or  the  dimension  of  the  boundary  zone,  can  be 

-3 

obtained  by  equating  the  corresponding  term  to  e  ,  which  is  approximately 
0.05  in  value.  The  lower  of  the  two  roots  m^,  is  given  by 

_  (n\-  ft* 73™"')*  i  os) 

"l  *  V  2 u  '  c 


for  orthotropic  materials.  The  calculation  of  decay  length  is  illustrated 
below. 


>  79) 


x.  is  decay  length  and  is  obtained  as  1.682  H.  For  isotropic  materials,  an 
a 

approach  similar  to  the  above  is  followed  starting  from  Equation  (72).  x^ 
is  obtained  as  0.410  H.  It  is  independent  of  material  constants. 

Equation  (79)  indicates  a  strong  dependence  of  the  decay  length  on 
the  material  properties.  The  decay  length  defines  limits  for  application 
of  a  decay  type  solution  of  the  form  given  in  Equation  (76).  A  minimum  of 
two  decay  lengths  is  required  for  the  beam  to  be  considered  long  so  that 
corrections  of  the  type  in  Equation  (77)  are  applicable.  Beams  with 
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lengths  less  than  2x^  must  be  considered  short  and  a  full  solution  of  the 
type  in  Equation  (71)  should  be  used  in  the  computations  of  refined 
clamping  model  solution. 

An  estimate  of  the  clamped  end  warping  cs placement  may  be  obtained 
by  integration  of  the  constitutive  relation 


u 
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u  is  nondimensionalized  with  respect  to  maximum  deflection  at  the  center 
according  to  Bernoulli-Euler  theory  for  convenience  and  is  expressed  as 
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u  0  is  the  axial  displacement  at  the  clamped  end  due  to  C2  restraint. 

V4 
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max  B-E 
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W  _  is  the  maximum  deflection  at  the  center  according  to  Bernoulli- 
max  B-E 

Euler  theory  and  is  given  by 


W 


max  B-E 
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(83) 
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The  clamped  end  transverse  displacement  w  is  of  secondary 
importance  for  the  present  study.  It  is  usually  Poisson's  effect  and 
therefore  much  smaller.  It  may  be  obtained,  in  a  similar  manner,  by 
integrating 


The  result  is 
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The  clamped  end  displacement  components  are  shown  in  Table  1.  It 


can  be  concluded  that  the  approach  adopted  in  the  present  study  provides  a 
simple,  reliable  model  for  ideal  clamping. 


Since  the  stresses  do  not  satisfy  compatibility  relations  exactly, 
the  above  approach  is  not  a  unique  way  of  determining  the  displacement 
components.  It  is,  however,  rational  and  the  most  direct. 

Concluding  Remarks 

On  the  basis  of  the  results  presented  above  the  following 
conclusions  are  reached. 

1.  A  way  of  determining  the  boundary  zone  stresses  is  presented 
which  indicates  that  the  interior  solution  is  best  represented  by  C2  model. 

2.  It  is  demonstrated  with  confirmatory  results  that  C4  model 
clamping  is  extremely  good. 

3.  An  estimate  of  end  zone  correction  for  long  beams  is  provided 
through  LDP.  The  region  of  local  effects  is  quantified  with  the  aid  of  the 
decay  length;  it  is  shown  to  be  a  strong  function  of  material  properties. 


CHAPTER  VII 


VALIDATION  OF  THE  THEORY 

Preliminary  Remarks 

It  is  demonstrated  in  Chapter  V  that  the  new  equations  yield 
exact  results  or  results  that  are  indistinguishable  from  exact  for  the 
static  examples  considered.  The  loading  is  uniform  in  all  but  one  case.  A 
linearly  varying  load  is  considered  in  the  exceptional  case.  The  validity 
of  the  theory  for  arbitrarily  varying  load  remains  to  be  established.  This 
is  accomplished  by  completing  the  following  three  tasks: 

(1) .  A  thorough  analysis  to  determine  error  estimates  for  the 

equations  is  presented. 

(2) .  Reissner  plate  equations,  specialized  for  planar  bending,  are 
demonstrated  to  be  obtainable  from  the  present  equations. 

(3) .  A  quantitative  demonstration  for  a  classic  benchmark  problem 
is  provided. 

The  benchmark  problem  is  the  response  of  a  simply  supported  beam  tc 
a  sinusoidally  distributed  loading.  Exact  solutions  for  this  problem 
appear  in  References  52  and  53,  which  facilitate  a  critical  comparison. 
This  is  a  generic  problem  which  has  been  used  as  a  test  case  by  others. 
Predictions  of  the  present  theory  are  compared  with  the  exact  solutions 
using  an  approach  which  yields  the  range  of  validity  of  the  theory  as  a 
function  of  beam  length-to-depth  ratio.  Consequently,  a  direct  indication 
of  the  applicability  of  the  theory  in  this  nonuniform  loading  situation  is 
obtained  for  specific  geometrical  and  stiffness  characteristics. 


The  stresses  given  by  the  present  theory  are  approximate  for 

nonuniformly  distributed  loadings.  It  is  desirable  for  the  errors  in  the 

equations  of  equilibrium,  compatibility  equations  and  displacements  to  be 

consistent  with  the  level  of  approximation  of  the  stresses.  To  study  this 

issue,  a  systematic  order  of  magnitude  analysis  has  been  undertaken.  The 

approach  adopted  and  the  underlying  philosophy  of  the  arguments  presented 
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are  similar  to  those  employed  by  Koiter  ’  in  conjunction  with  a 

critical  study  of  shell  theory  equations. 

The  magnitudes  of  spatial  derivatives  are  estimated  in  the 
following  way:  a  wavelength  for  the  deformation  is  defined  such  that 


is  the  maximum  absolute  value  of  the  quantity  <t>  in  the  region  under 
consideration.  X  is  associated  with  the  wavelength  of  load  variation  and 
the  wavelength  of  deformation.  Derivatives  with  respect  to  z  are  dealt  with 
in  a  similar  way. 

This  implies  that  the  smallest  wavelengrh  of  deformation  to  be  considered 
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Let  L  be  the  measure  of  beam  length.  For  applications  of  interest 
here,  L  =  0(X)  and  H/X  is  small. 

In  bending,  0 ^  is  the  largest  stress  and  is  chosen  as  a  convenient 


reference.  Let  0  be  its  maximum  value.  With  the  aid 

of  Equations  (22), 

(16)-(18),  the  following  estimates  are  obtained: 

1 

1m1  “  0 (o  H2) 

(87) 

|q  |  -  0 (o  H2/X) 

(88) 

• 

|0  |  -  o(^) 

H 

(89) 

Equations  (23),  (24)  and  (88)  permit  the  estimation  °faX2  and°z2* 

|<J  I  ~  0(0  H/X) 
xz 

(90) 

1 

|o  j  -  o(  <jh2/X2) 

(91) 

•  zz  ' 

3 

The  nonclassical  axial  stress  is  represented  by  the 

(33A).  From  the  result  (88),  this  may  be  estimated  as 

underlined  term  in 

. 

°<kx0BV> 

(92) 

The  subscript  "NC"  refers  to  the  nonclassical  part  of  the  stress. 

Error  Estimates  for  the  Equilibrium  Equations 

In  the  classical  theory,  the  stress  equilibrium  equations  are 
satisfied  identically.  Due  to  the  nonclassical  bending  3tress  in  the 
present  theory,  stress  equilibrium  is  not  satisfied  exactly  for  nonuniform 
loading.  In  order  to  facilitate  the  argument  that  follows,  the  stresses 
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0  and  0  are  written  as 
XX  xz 


XX 


<7  -a 

xz 


+  0  „„ 

(93) 

xxC  xxNC 

xzC  +  °xzNC 

(94) 

The  subscript  "C"  refers  to  the  classical  part  of  stresses.  It  has  been 
14 

shown  by  Seewald  that  0  contains  higher  order  terms  in  addition  to 

xxNC 

the  term  used  in  the  present  theory  for  nonuniform  loading.  The 

H2 

representation  for  0  is,  therefore,  valid  up  to  terms  of  0(0— r-k  ).  The 

XX  X 

nonclassical  term  is  seen  to  be  a  function  of  beam  geometry  and  the 

relative  stiffness  represented  by  k^.  Equation  (33A)  for  0  provides  a 

2  2 

good  approximation  for  small  values  of  the  combination  H  A  .  The 
nonclassical  axial  stress  effects  are  more  significant  for  orthotropic 
materials  with  large  kx  values. 

Introduction  of  Equations  (93)  and  (94)  into  (25)  leads  to 


axzNC,z  +  Q,xx  ~I  (z3  '  I  c2z) 


(95) 


The  above  permits  to  estimate  the  nonclassical  shear  stress  to  be  of 

3  3  H2 

0(k o H  / X  )  and  the  error  in  equilibrium  is  of  0(0  — r  k  ).  It  can  be 
x  x3  x 

concluded,  therefore,  that  the  error  in  the  stress  field  is  at  most  of 

H3 

0(j  — r  k  ).  This  is  consistent  with  the  original  approximation  for  0  . 

X3  x 

Error  Estimates  in  the  Compatibility  Equations 

The  relevant  corapatibil ity  equation  expressed  in  strains  for  planar 


binding  is  given  by 
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e  +  e  -  Y  =0 

XX. zz  zz.xx  xz,xz 


It  is  convenient  to  express  Equation  (96)  in  terms  of  stresses. 
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Substitution  of  Equations  (33A),  (94),  (23)  and  (24)  into  (96A)  leads 
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The  error  represented  by  the  underlined  term  is  of 
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the  nonclassical  shear  stress  is  of  0|O  — x - - - 

V  X  '33  <4  -  «„> 


It  is  concluded,  therefore,  that  the  error  in  the  stress  field  is  at  most 


of  0  |  a  —3 


This  is  consistent  with  the 
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original  approximation  for  0^  in  Equation  (33 A)  and  with  the  error  estimate 
found  for  the  equilibrium  equations  based  upon  Equation  (95). 


Error  Estimates  in  Displacements 


The  error  in  u  due  to  the  error  in  stresses  may  be  estimated  from 


equations  (7)  and  (27), 


u 

* 
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(27A) 


The  error  in  W  may  be  estimated  from  Equations  (7)  and  (29) 


w  =  -  u 
,x  ,z 


(29A) 


H 

The  error  in  the  stresses  is  of  0(<7  — r  k  ). 

X3  x 

(27A)  and  (29A) 


Consequently,  from  Equations 
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In  view  of  Equations  (98)  and  (99)*  it  is  justified  to  omit  terms  of  the 

order  indicated  or  higher  in  the  expressions  for  displacements  in  Equations 

H3 

(31)  and  (32).  The  underlined  terras  in  these  equations  are  of  0(- —  —?) 

0  H4  .  11  X 

and  0(- —  — -),  respectively.  It  is  consistent  to  ignore  these  terms  on 

11  X 

the  basis  of  the  above  discussion,  so  the  approximations  made  are 
consistent . 

Summary 

The  study  of  the  order  of  magnitude  of  the  errors  in  stresses, 
displacements,  equilibrium  equations  and  compatibility  equations  has 
demonstrated  that  the  present  equations  are  self-consistent  and  provides  a 
valid  approximation  when  the  error  terras  are  negligible. 

This  implies  (H/X)  is  sufficiently  small. 
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Relation  to  Reissner  Theory 

Reissner's  plate  theory  equations  are  derived  by  using  a 

37  38 

complementary  energy  principle  ’  .  The  definitions  of  the  kinematic 

variables  are  clarified  in  the  latter  paper.  Weighted  kinematic  variables 
naturally  arise  due  to  the  approach  used  to  develop  the  theory.  For  planar 
bending,  they  are  given  by 
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(101) 
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In  the  following  derivation,  the  Reissner  variables  are  constructed 
using  present  theory  displacement  expressions.  Equations  (100)  and  (31) 
permit  w  to  be  written  as 


W  »  W  - 


3V13M 
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<f>  is  identical  to  $  ^  £iven  Equation  (37). 
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Equations  (104),  (37A)  and  (18)  allow  the  following  to  be  written. 


♦  ■ -  5  * 


(103) 
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The  above  is  the  Reissner  relation  for  transverse  shear  strain.  Intro- 
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duction  of  Equation  (102)  into  (37A)  and  the  use  of  Equation  (18)  results 
i  n 
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(104) 


By  virtue  of  Equation  (103),  (104)  may  be  rewritten  in  the  following 

familiar  form: 


13 
A  ' 
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The  above  demonstrates  that  Reissner  equations  can  be  obtained 
from  the  present  theory,  a  fact  which  further  establishes  the  validity  of 
the  new  equations.  Also,  it  can  be  observed  that  all  of  the  essential 
physical  effects  do  not  appear  in  the  Reissner  theory.  This  is  because  a 
knowledge  of  Reissner 's  variables  does  not  permit  the  determination  of  the 
response  throughout  the  structure.  The  effect  of  nonclassical  axial  stress 
on  response  is  totally  lost  in  the  averaging  process.  However,  by  the  use 
of  the  relations  presented  above,  the  response  in  terms  of  Reissner 
variables  can  be  converted  to  obtain  the  response  throught  the  structure 


Beam  Under  Sinusoidal  Loading 

The  problem  under  consideration  is  described  in  Figure  10.  The  two 
dimensional  elasticity  solution  for  it  is  given  in  Reference  52  for  an 
isotropic  material  and  in  Reference  53  for  an  orthotropic  material.  A 
solution  to  the  above  problem  has  been  obtained  by  using  the  present 
theory. 


One  half  wavelength  of  the  deformed  beam  is  isolated  for 
consideration  and  is  treated  as  being  simply  supported.  Although  only  one 
half  wave  length  is  considered  the  results  are  applicable  to  the  case  of 
general  loading  of  the  form  sin  This  follows  from  the  fact  that  each 
half  wavelength  may  be  considered  separately  with  an  appropriate  reduction 
in  beam  length.  The  coordinate  axes  and  notation  are  also  given  in  Figure 


The  boundary  conditions  to  be  enforced  are  given  in  Equations  (39). 
In  addition,  the  following  are  also  satisfied: 


O  (0,z)  *  a  (L,z)  -  0 

XX  XX 


(106) 


The  transverse  displacement  component  W  for  the  above  boundary 
conditions  is  obtained  by  integrating  Equation  35.  The  result  is 
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The  stresses  O  ,  0  and  o  at  any  section  x,  are  given  by 
XX  xz  zz 
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It  is  possible  to  obtain  further  refinements  to  the  present  theory 
by  retaining  the  underlined  terms  in  Equations  (31)  and  (32).  The 
constitutive  relations  are  obtained  as 
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The  axial  stress  distribution  is  given  by 
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The  corresponding  transverse  deflection  is  given  by 
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The  underlined  terms  in  Equations  (111)  and  (112)  are  the  refinements  to 
the  present  theory.  Equations  (109)-(110)  will  be  called  Approximation  1 
and  Equations  (111)  and  (112)  are  Approximation  II  in  the  following 
discussion.  The  purpose  of  the  above  refinements  is  to  determine  the 
influence  of  the  underlined  terms  in  Equations  (31)  and  (32)  for  this 


i 


74 


problem. 

Results  and  Discussion 

The  results  are  presented  in  a  common  format  and  appear  in  Figures 
11-18.  For  the  orthotropic  beamf  the  following  properties  are  chosen: 
E^  =  25,  =  =  0.5  and  «*  0.25.  These  are  the  material 

constants  used  in  computing  the  exact  solution  in  Reference  53. 

The  relative  merits  of  each  theory  under  consideration  are  assessed 
on  the  basis  of  the  percentage  error  with  respect  to  the  exact  solution. 
For  the  present  purposes,  a  five  percent  error  is  assumed  to  be  an 
acceptable  limit.  The  range  of  beam  length-to-depth  ratio  in  which  the 
error  is  less  than  five  percent  is  considered  the  range  of  validity  for  the 
theory.  The  point  at  which  a  theory  just  exceeds  the  limit  is  a  cut  off  or 
limit  value  of  beam  length-to-depth  ratio. 

The  salient  features  of  the  results  are  listed  below. 

(1)  Present  approximations  provide  superior  predictions  for  the 
response . 

(2)  The  improvements  are  more  significant  for  orthotropic  beams. 

(3)  The  Approximation  II  appears  to  be  only  marginally  better  than 
Approximation  I.  Approximation  I  is  fully  adequate  for  most  applications, 
there  fore . 

In  Figure  11  a  curve  that  corresponds  to  Reissner  theory  equations 
(103)  and  (105)  is  also  shown  for  comparison.  The  Reissner  theory 
prediction  appears  to  be  in  excellent  agreement  with  the  exact  solution. 
However,  this  is  illusory  as  the  quantities  under  comparison  are  not  the 
same.  Reissner's  displacement  variable  is  a  weighted  average. 


PERCENTAGE  ERROR  IN 
MAXIMUM  TRANSVERSE  DEFLECTION 


Figure  11.  Percentage  Error  in  Maximum  Transverse  Deflection  for  a 
Sinusoidally  Loaded  Isotropic  Beam 


PERCENTAGE  ERROR  IN 
MAXIMUM  TRANSVERSE  DEFLECTION 


PERCENTAGE  ERROR  IN 
MAXIMUM  AXIAL  TENSILE  STRESS 


PERCENTAGE  ERROR  IN 
MAXIMUM  AXIAL  COMPRESSIVE  STRESS 
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Figure  14.  Percentage  Error  in  Maximum  Axial  Compressive  Stress  for 
a  Sinusoidally  Loaded  Simply  Supported  Isotropic  Beam 


Figure  15.  Percentage  Error  in  Maximum  Axial  Tensile  Stress  for  a 
Sinusoidally  Loaded  Simply  Supported  Orthotropic  Beam, 

®  1 1  /G i ^  “  50,  m  ^ 
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PERCENTAGE  ERROR  IN 
MAXIMUM  AXIAL  COMPRESSIVE  STRESS 
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Figure*  17  and  Id  chow  the  axial  stress  distribution  through  the 
depth  of  bee*  at  midspen.  The  agree* ant  of  Approximation  II  with  the  exact 
ia  excellent  uhile  Approximation  I  can  be  considered  satisfactory.  The 
effect  of  neoclassical  axial  at res*  ia  seen  to  be  more  pronounced  for  the 
orthotropic  beam.  Tensile  stresses  are  predicted  by  both  Approx  mat  ion  I 
and  Approximation  XI  near  the  center  of  the  cross  section  as  seen  in  Figure 
19.  There  ere  errors  in  the  present  spproximat ion  in  this  portion  of  the 
cross  section.  Maximum  stresses  are  predicted  quite  well. 

It  is  interesting  to  note  that  a  Re  is soar  or  Timoshenko  type  theory 
would  have  given  the  same  result  as  classical  theory.  Monctassical  axial 
stress  contribution  is  not  present  in  those  theories. 

The  transverse  shear  stress  and  transverse  normal  stress 
distributions  at  sections  where  they  are  maxi  mm  are  presented  in  Tables  2 
and  3.  The  agreement  with  the  exact  values  is  satisfactory. 

Pome I us tons 

A  qualitative  validation  Is  provided  through  a  consistency 
analysis.  The  theory  is  farther  established  by  showing  that  1* issuer 
theory  of  plates,  reduced  to  planer  beading,  can  be  developed  fro*  the 
present  theory.  A  quantitative  validation  is  provided  Uirovpi  a 
correlative  study  with  the  exact  solution  to  n  classic  benchmark  problem  * 
-  tbe  response  of  e  siapty  supported  beam  to  e  sinusoidally  distributed 
loading. 
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The  abow  cannot  b«  Mtilfiad  identically  for  all  a.  In  view  of  the 
appro  stoat  ion  inherent  in  the  stresses,  it  ia  desirable  to  satisfy  Equation 
(127)  to  the  same  degree  of  approximation,  the  terms  may  be  rearranged  in 
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»y  forcing  the  first  term  in  brackets  to  sero,  the  following  equation 

of  motion  ia  obtained 


M  -  Q  -  D  It.  (115A) 

t*  * 

Where  #,  ia  the  rotation  related  variable  given  in  Equation  (36).  Since 
*  -2 

the  «esder lined  terms  are  at  most  of  0(  0  ~  k  ),  Equation  ( 1 15A)  implies 

* 

that  the  approximation  inherent  in  the  equation  of  motion  is  consistent 

with  that  of  the  stresses  and  compatibility  equation.  This  is  similar 

to  the  static  siteatioa  of  Chapter  VII.  Equation  (11 5A)  is  the  overall 

23 

equation  of  motion  of  the  original  Timoshenko  theory  .  Equations  (23), 
(117)  and  (lit)  field  the  secood  equation  of  motion  ss 
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The  eaiet  itrtit  fitltibatioo  it  ottoioo*  vitk  the  oi4  of  tguetiooo 
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As  («•  ,  Hoot  ion  ( 140)  |iwi  m»  lor  do  phene  vvlac  ity 

poroneter  (.  Dm  tooor  roo<  wrtip-li  to  floaorol  ww  velocity  lor 
to  rose  ic  «•«».  Dm  looor  knack  Is  plot  i *4  age last  Dm  amU^tk 
percenter,  (  lo  flforo  It,  (or  onset,  pro  root  ao4  Ttnkwkc^  dootiw, 
Poisson's  rot  lo  lo  tohon  oo  0.29  lo  Dm  oonpetet  lone .  TW  oooct  oetotloo 
koo  Woo  rocoootrectod  Iron  Dm  analysis  to  toforooco  H.  Dm  nm  for 
Tlnoshenko  theory  lo  kitil  apes  Dm  original  aDoUoaoU,aot  Dm  lot  or 
versions  (Aid  "ad)eec”  tW  sheer  oorroccioo  for  tor .  Dm  oooct  eolation  to 
(or  oo  infinitely  old*  plot#.  Dllo  corresponds  to  ptooc  otrolo.  C  oM  v 
ore  therefore,  replaced  by  t/(  I  -  v  *)  and  v/Cl  *  *)  lo  DM  pro  root 
•yoiioM.  teoelto  toe  to  Dm  present  theory  ore  lo  eocelloot  at  rv  tenet 
vttk  Dm  eooct  oetotloo  for  Dm  roots  of  aeve lengths  considered.  Dm 
ooya^totic  eoteeo  of  DM  phone  velocity  poroneter  ore  0.9 290,  0.H12  ooO 
O.(ll)  for  eooct,  pro  root  and  Ttoeehenfco  theories,  respectively. 

DM  oooct  eolation  loots  to  infinite  hr one he s  (or  the  dispersion 
carve.  All  the  hither  bronchos  tench  ooyptoticolly  the  cheer  neve 
velocity,  fenf insert 0|  bend 1st  theories  glee  only  too  brooches,  the  second 
of  vhich  is  o  poor  oppron toot  ion.  Dm  second  osynptote  for  the  phase 
velocity  poroneter  is  1.470  end  l.Af  for  present  end  Ttoesheoho  theories. 

Dm  ooynp  totes  for  the  phone  velocity  poroneter  ohtoined  as  tot 


c/c. 


tiepPoa  mft  Uti««  (kMT|  «•  0.9411  m4  1.4419  far  tli  (ini  m4  nwK 
mull  rMpNtivlf. 

Ita  rmlu  cm  Pa  «m4  m  (artltr  aelidete  di  cm  hmuwi . 

ft*  *»mrt  i«  aealapeaa  l*  ft*  «(ft*ti*a  Nftf  *f  ft*  tutu  ft*ff.  Ac 
nail  aalatim  i*  m  artPecraptc  *l*k  t*  ***M  far  tin  parpaaa.  It  U 
ft**l*H  Ic  Appmdic  A.  leatrapic  rmlu  crc  ftuiw*  Of  apaciel  lutiw. 

I***ltt  cl  ft*  «ll4(i*i  acadp  fcr  m  laacrapic  *lc0  ipyiir  Ic 
lifM**  1ft- 22.  Hm  parcactapa  crvcrc  ic  pPaaa  cclcclif  aPam  Ic  li|nr*  20 
tadlaeta  iIm  ftc  pwim  tic  cry  predict tarn  ere  cmpcrcOlc  u  ItepOm  aft 
Leri  cccc  ftiirf  reed  cat  tail  ere  cal  1ft  Ic  ft*  act  ire  raape  cf  etaaOeracM 
pircaw  earn i dared.  fl>aii  tt  aeft  22  iadimee  ft*t  ft*  praam*  toaary 
Oaedlap  ecreee  aaclmtaa  arc  amerler  tc  Permelti -taler  aft  TfcmtPmOa 
cOeared  cat  read  cat  fta  add  acreaa  die  triPet  im  ftmft  eba  depeP  la  ic 
acacllaM  asnaMM  alft  fta  aaacc  actcaias.  Sceyfcm  aeft  tedmm  tliary, 
Ic  lea  praam*  fare,  amme  predde  ami  id*  acreaa  aaclmtaa. 

Oarraapacftlat  rmlu  far  m  erCPecrepla  atcO  era  prat cat ad  Ic 
lipftM  21-25.  tidier  PaPederd  tree**  ca  ft*  laacrapic  ecaa  era 
■Octree  A.  freeeat  lOaary  predict  iem,  Pemcet*  aim  mra  pr  aaacc  eat 
iftpraeaamca  am  tOaaa  at  TtmePeaPe  aeft  it  race!  It -talar  ftcaeiat.  far 
caaaple,  Occaft  «pm  tOe  par  tea*  ape  arrar  ic  tPa  add  acreaa  aOaaa  la 
flpera  24*  Tima  Paata  ftaary  ia  taliA  ap  *#  L/V  •  15.5  aoeraaa  fta  praam* 
ePaary  ia  a  peed  appradmeim  ap  *a  l/l  •  2.5.  Praam*  tPecry  pradeCim 
af  tPa  adal  atram  AiatrfPatim  tPteagl  AaptP  ia  ia  acealtm*  apraamat 
dlP  tPa  acme  ratal c.  ftp  caairaat,  aa  aParn  ia  fipar*  25*  TimaPrefcc  act 
tataoalli -talar  praftietiam,  aPicP  art  Umar  tlrwpi  ft*  AaptP*  arc  ia 
eary  p car  apraamat  dtP  *Pe  cam*  mlatim. 
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No  rooolto  for  StopMo 


Uoi 


iMorp  oro  •«  if  M« 


km  fmloH  ooly  f«r  iootropic  Mirfiil*. 
likutiw  hokloi 

Vikrrtin  Mkooior  for  kom  mill  ooriooo  o»0  rootroioto  »• 
okuioaf  ooiog  tkr  procotforo  ootliood  io  Mforooooo  >>  mi  SI. 

Typical  rooolto  for  o  ctaapod  orfknrofu  km  appoor  i«  fi|*r#«  H 
aod  It.  «*  (  roforo  to  (No  ftorooolt t -Color  frogooocy  for  «  tioflr  ooo^oriol 
km.  Tkt  fifftroMM  moot  *No  cooirot  ooomt  rifit  ptM«  or# 

■or*  (Mo  (Moo  rtflKiH  io  iM  frMoaaty  rolio  grapfco. 

1M  rooolto  for  iootropic  kom  ami  for  itaplf  oopportod  ood  cl«H 
rootroioto  otoo  Moo  Moo  ottoiood.  TM  rooolto  Moo  oMm  |mt*l 
agroomnt  of  troodo  oitk  TtaooMoM  (Mory.  TMro  io  ooooitiotty  *• 
aodoliog  of  (M  N ooo dory  rootroiot,  Mt  It  to  loot  prooooocc 4  (Mo  *‘*i 
(mf  io  tM  ototic  opplicatiooo. 
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Control  Meant  Ratio  for  a  ClanpeR  OrtRotrojic  Re 
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ckthi  u 

AM  SLMMDTTAXT  tOCKUMC  TMSOtT 
fwlUt-ry  tsosrts 

la  this  dMpttr(  so  ilnaurjr  hwklii|  theory  i«  4mleH  with  ch* 
•i d  of  the  (fatioM  fmloH  ia  Chapter  IT.  Tho  primary  obJ«eti*t  of  the 
fraaaat  ataiy  to  to  cooper*  beck l lot  load  prodietioes  with  chooo  of 
Tteoshsoko  aod  Boroool  11 -Color  Choorioo.  Coaseqesotly,  linearised 
•fOtioM  or*  etillaed.  A  folly  ooeliaesr,  lorft  disploceasnt  theory  l« 
bsyaod  tho  seep*  of  tho  present  writ. 

At  tho  oooot  of  beck  Hof,  Cqestioa  (22)  io  «N«wd  to  protlft  first 
approolaatico  to  tho  axial  otrooo.  Tho  roasioiot  stresses,  which  oro 
cooolotoot  with  thio  spproc  last  too,  oro  taboo  to  ho 

"  «  2  2 

O  n  •  “J|  Cc  *  l  )  ( 23A) 

»„  *  iF‘7  <*»> 

Tho  fora  for  «M  io  slightly  different  froa  that  of  the  correspond i a* 
static  distribotion  hocaoao  of  difference*  io  the  hoaodary  condition*. 
Stroaa  fro*  cooditiooa  oo  c  •  ♦  c  aorfaco*  aro  satisfied,  to  this  respect, 
it  ia  aoolofooa  to  tho  dyooaic  stress  distrlbatioa  fivsa  io  tqastioa  (24A) . 
Also,  tho  strosaos  aro  copra  seed  io  toras  of  tbs  hoodioc  oo— t  M  and  its 
derivative*  ooly.  This  is  hocaoao  fqestioa  (1C)  is  no  looter  valid,  hot 
tho  special  distributive  of  stress  over  tho  cross  soctioo  roosins  sioilsr 


life 


to  clko  cUmIu!  com. 

•jr  following  tbs  dtwloyMit  bm4  is  Qupttr  IV*  tin  t^rtNioa  for 
tho  bonding  mommt  U  obtained  u 


b^  It  •  for— tor  defined  io  Mention 


s 


(142) 


(US) 


(142) 


lbn  buckling  is wot ion 

An  oooroil  equilibria*  •fotln  U  derived  fro*  static  squill  brio* 
of  a  4s  fisc  tod  boa*  olso*t  according  to  Cbo  ad  jacsot  equilibria*  approach. 
fipiN  2t  ehowe  tho  for coo  oa  a  do floe cod  boa*  oloaoat  of  fiaito  length. 
Tbs  *o*aat  equilibria  require* 

N  •  PW  ♦  N  (144) 

o 

M  io  cbo  coaocoac  aoaooc  at  a  •  0*  abicb  io  oood  to  satisfy  particular 
o 

boundary  cooditioaa.  With  tbs  aid  of  tquatioaa  (143)  aad  (144),  cbo 
following  equation  for  the  bwcklod  ooofigu ration  is  obtainodt 


« 


•h1  u  *i  S77T 


> 


0 


(145) 


tquotioa  (143)  io  soon  to  bo  of  the  oa*o  for*  as  tho  classical  colu*n 
equation,  which  *ey  bo  obtaiaad  by  sotting  a#  •  a(  •  •  0  in  tho 

a bows. 


Por  eiaply  aapportod  oada,  tbs  boundary  conditions  to  bo  ooforcod 


» 


III 


ere  given  is  Bqaat loss  Of).  A  bock l las  ratio,  P^Pg,  is  obtained 

with  the  aid  of  Meatless  Of)  aad  (145). 


vr«  - - ~-B- 

'*  ’  f** 


(144) 


fQ  is  chs  beck  Hag  load 

eiaply  nfyortH  col ana. 


Fj  is  dM  clsssicsl  Calsr  heckling  load  for  s 


•  s,  .1 


*«  " 


(147) 


p  •  V  I/A  is  the  radios  of  gyration  of  tks  boss  cross  ssctlas.  Tbs  ratio 


for  s  definition  of 


L/p  is  dM  eleadnraasa  ratio  of  ths  eolsas. 

B< «■  ties  (144)  offers  a  convenient 
s leads rasas  of  a  col  ran. 


Pif  -  - -f  -j 

^  1  1  ♦  tVaJ 


( 1474) 


s  Is  a  universal  stands rases  par vaster  da  flood  by 


•WT 


(144) 


This  psraaacsr  Is  s  faoctioa  of  both  gaonatry  sad  etiffaess  of  tbs  col  ran 
asterisl.  As  its  vales  increases,  clsssicsl  Beter  coluso  behavior  is 
approached .  A  "short"  colrao,  thea,  inpliss  s  deparmre  fros  ctsssical 
Baler  behavior  and  is  characterised  by  suit  valves  of  s. 
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toulti  ad  tiic»Mio» 

Th#  beck 1 lag  load  ratio*  art  calculated  for  aa  ortkotro^ic  tiapl/ 
e up ported  co loan  with  I^/C^  •  30  aad  v(J  •  0.).  fof*g  prediction* 
according  to  do  pnmt  aad  Tlaoaheako  theory  are  0.529  ao4  0  429 
reepectively,  for  l/a*  20.  Tlaoaheako  theory  five*  a  lower  value  for  the 
buck 1 log  load  ratio  aad  la,  therefore,  overly  cooaervativ*.  Departure* 
fro*  claaeical  theory  are  aehataatia)  ia  thia  cave. 

If  a  five  percent  depart  ore  froa  claaeical  theory  to  aet  aa  a 
practical  liait  to  differentiate  betooea  a leader  aa d  abort  col  wane ,  thee  a 
threahold  vale*  a^  for  the  miverael  a lead* me aa  paraaeter  cae  be  obtained 
froa  tguation  (I47A).  A  aaefal  approaiaatioa  ia 

*t  •  15  C 149) 

The  correapoadiag  poaetrUal  alaadaraaaa  ratio,  L /p  ,  value*  for  aa 
ieo tropic  aaterial  with  v  •  0.)  aad  aa  orthotropic  aaterial  with 
«,l/Ci3  -  10  -d  Vjj  •  0.3  are  20.0  aad  90,  reapectively.  Thee*  value* 
ladicate  that  geoaetr ic  ateaderaeaa  atooe  ia  aet  iadicatlv*  of  coluan 
behavior. 


II) 


aunta  x 

H M  pn)lw  mmcU(M  ate*  Lytrotlw—l  offocto  is  fibtr 
rat  of  or  cod  mil  MCrU  cm^mUm  hra  b«M  MrUtr.  Om  of  ebo 

«r<  aarioao  ftf—m  of  hygncfemil  eooditioaiot  id  do  dogrodatioo 
of  otiffoooo-volotod  a4  om|rt  tolotod  proper t loo.  to  dw  ywiwt 
otady,  Nfloiif  !o  ploool  M  pro party  dagradotiaa  offocto  U  otiffaaoo 
critical  oppliootiooo.  taolllot  otrooooa  Am  «  troooiaot  Mlitm 
fitirikadoi  oro  oat  otadiod. 

Sot  Loot oo  of  loot  of  parforooaoa  Am  to  hytrocterMt  oooditioaiag 
for  aoidiraactaool  eaopoaitaa  oro  pro— tod  lo  fblo  cOapcar.  Ortbotropic 
oolotiooo  dorioad  to  proriooo  dwpcoro  oro  otiticod  to  ofccoio  ckooo 
oociaotoo.  Dm  fotloolog  oolotiooo  oro  applicabla  to  Midiractioaol 
loyapa.  Dm  ootioocoo  —aid  Lo  rloood  oo  proridios  —  litotioa 
ioforaot  too  for  taoigit  ioto  pot  oat  to  1  practical  coaoo— acoo .  It  io 
ooticipotod  tlot  ooopootto  otroctoroo  otU  off  ooto  ply  loywpo  will  adhibit 
groat or  bygro  thermal  offocto. 

a 

Hod-tool  Eroportioo 

Tbraa  oaoiranooatol  eoaditioao  liotod  be  loo  oro  oo  lac  tod  for  tht 


otaiyt 
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CooGittoo  At  Dry  300*? 

CoMitua  Dt  1.0*1  Viator#  Cnimi  300*? 

CoMitioa  Ct  1.41  Viator#  Cm(M(  300*? 

CooGii i«  a  it  (or  ftftriK*.  CoaiiiiM  •  rapratoo tt  <  rvilittic 
•iiMiiM.  TV  t^ial  l«t*l  of  viatara,  1.0*1  4?  w»i#tc.  it  • 
appro aioat  too  to  loop- taro  aircraft  aaratea.  CooGtctoa  C  ta  aatoractoo  for 
ao  At  MON  Graph itt/Bpoay  aoG  rapraaaoca  cha  ooat  aoaora  OapraCat too 
taook  at  tha  aalactaA  rooparatora.  Typical  proportiaa  have  haao  choaoo 
ait 4  tha  ait  of  aaparimotal  Goto  proviGaG  io  Bafaraoca  *f.  Thoaa  era 
praaaotaG  io  Tohla  4. 

Tha  raapooaa  of  ao  ortho tropic  koto  aoGar  alaaaoiGally  GiatrihataG 
looGiag  ia  coapatoG  at  tha  thraa  tavata  of  hygrothoroal  cooGitiooiap.  TV 
raaalta  ara  praaaotaG  ia  Piparaa  *1  ood  32.  Tha  raapooaa  oGcr  eooGitiooa 
A  aoG  •  ia  a  littla  Gifforaot,  Vita  cooGicioo  C  raapooaa  aahihita 
|  aobataotial  Glfforaocaa  froa  A  ooG  t.  Tha  axial  atraaa  Giatrihatioo  at 

tOGapao  thoaa  io  Fipara  32  laGieataa  aceaotoatat  oooclaaaieal  affacta 
aoGar  hyprotharoal  cooGitiaoiop.  fey  oeatraac ,  tha  ctaaaicol  or  a 
I  Thaoahaaho-typa  thaory  faila  to  praGict  hyprotharoal  aaoaitirity  of 

haoGiop  atraaa.  Tha  taoaila  atraaaaa  proGictaG  oaar  the  contra  of  ct-oaa 
aactiao,  ahoaa  tha  oaatral  uia  arc  Gao  to  orrora  io  cha  prcacot 
I  appro xioat ioa  ia  thia  portioo  of  tha  croaa  aoctioo. 

ClaopcG  orthotropic  boao  raaalta  aheoiop  the  frequency  ratio  chanpea 
Gao  to  hyprotharoal  cooGitiooiap  appear  ia  Pigwre  33.  IV  claopiap 
*  condition  corraapooGa  to  C2  rcatraiot.  TV  behavioral  treoGa  arc  aioilar 
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Ta4l«  4.  NtfiriiM  •(  dM  ftrttxwyii  HMru!  »i  OnliiMM  &,  1  m4  C. 
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ff« «(•  MfaUttflM,  •  flM  flfWl  Mi  •  M 
fwr  OMiltiM  •  Mi  OMiftiM  C, 


rciMtlM  to  t 


»•  «Mti  •<  dM  At  am .  Caa<m«  C  Ami  mM(W  i«l 

•n«IM  l»wm  AffACCA  «  iMlUoi  Im4«  •(  *  Mlf) » 

mIw  iffNf  »•  fifir*  K.  TW  (mk«»h 

A  «M  •  Mt  W*  Atf*lft*AM  U  till  tWfU.  Cm4(Um  C,  Mr* 

(MMt  MflM*  NiffMM  WAultM*. 

in*  MMftt  i«|  tM4  nut  tm  A  KMtf  lAfUM  «dt«*  i<  <lM 

(MfalM  for  Mfi*M  )*•#)•  of  Millar*  «M*rn i«  M  um*iii *r* ,  TW 

•Mftul  mH  m  mi  n»i|«ii  i«  iiHmmi  •<  n*  «(f*i|Mi t«* 

IMfMItM  HI  MliMf«4  Will  (It  lU  Af  *1  MMlfftCal  MMi  *T*M«I*I  I* 

IfliNM*  |A  la  nil  M4*l<  a  |MliMl«  MrtMlMl  *r***rir  *1  n* 

■M  f 'is  Miff  lAl  AC  AAf  «M  d^MAItif  AM  MlllM*  W*lMl  it  •#**••*<  •« 
A  (MMliM  A f  MM  t <■»■*•« "«♦  *f  fMfATtf  AM  A  Mill  lift  ft  Af  (Aflff  ,  !%»*• 

I«l*f  <****«!  M  IM  *M  IMperMMf  PMIMIHI  A f  CAA  Mill***  AAAtAAt 
AM  IN  tlAAA  tf AAAiC tAA  IM^MAlM*.  TW  flMt  pMfffli**  AT*  AA««MM  t* 
f(Mi*  CM  MM  AM  CM  fTAAA  |*ApC*liM  *1  <AA»Olif  If*  f*M»<*<  If  «*i«| 
•ittMMAMilC. 

tM  tAAAlCA  AT*  pCA—ACM  i*  FifVt*  I).  TWf  iMicAtA  that 
IftMIIWfMl  AffACCA  M«A  AiMf  »  Afl  AAAAA*  «A  CM  lMlli«|  IaM  flti*  I* 
CM  tAAfA  At  l*MH*»*M«  10*f  *  I Vt*f  f A*  CM  ***»#*  At  Mi*t«r*  CAAtAAt 
CAAAi  A»*M.  At«AAC«ACM  AffACCA  CM  tA  Miltdf*  fit*  *f  At  A  »A*A  *aM  l* 
CM  ligMf  UA»*f  At  A*A  tAAfAA.  TM  Milt***  MlCtMt*  liAlt  i* 

in*tMCi«llf  MArH  At  AigMt  AAA  lAMpgf  AtAtAA. 

TM  fall**i*c  aaa  MMl*ti*M  W*aA  «a  tM  ilw*  *t at*: 

I.  Cm  At  Ml***  CttfkiltftfMf  *  Iff*  t  AAA  a*  ft*  AAtAtial*  in  tt*A 
aM*A  AfAtiCACiAAA  AO  AAt  f*»*  ***!***  ftfllCM  iA  AiAtol  At  AA  AitCtAft  *»<f( 


ill 

•itaftfi m  «ki(k  U  Hpwmu4  lyyinUy  Omnium  i. 

3.  9m  •!  f 1mm  utrUU  it  mmr  mtmmim  mmdiiim* 
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l»rt  29.  HyjTotherosl  Effects  on  the  Maxiaua  Deflection  for  a 
Siaply  Supported  Orthotropic  Bern  Under  Sinusoidally 
Distributed  Load inf 
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corclusims  in  RtooM*n*rta*s 

A  aw  engineering  dworj  for  planar  banding  haa  boon  developed, 
validated  and  appliad.  lea  predictive  capabilieiaa  have  baan  firmly 
aotabliabad  through  correlative  atudiaa  with  exact  oolutioaa  and  a 
systems tic  eona latency  analyeie.  The  theory  haa  boon  auccaaafully 
aodifiad  w  at  to  apply  to  dynawica  and  atatie  buckling.  Xt  accounta  for 
thraa  eaaential  phyaical  affacta  -  -  -  transverse  ahaar  atrain,  transverse 
normal  atrain  and  nonclaeaical  axial  atraaa.  The  equations  era  aa  aiapla 
to  apply  aa  any  tranavaraa  ahaar  da format ion- typo  theory  yet  they  provide 
response  throughout  the  structure. 

The  peasant  theory  predictions  are  superior  to  other  comparable 
engineering  theoretical  predictions.  The  theory  yields  exact  results  for 
the  case  of  uniformly  distributed  loading.  For  nonun i form  loading,  it  haa 
been  validated  by  means  of  a  thorough  consistency  analysis  of  a  qualitative 
nature  and  by  quantitative  correlative  studies  with  classic  benchmark 
problems . 

The  theory  includes  all  information  oontained  in  Reiaaner  plate 
theory  for  planar  bending  and  Stephen  and  Levinson's  dynamic  theory 
specialised  to  thin  rectangular  cross  section  beams. 

There  ia  a  sensitivity  of  the  predictions  to  boundary  restraint 
modeling.  It  is  more  pronounced  in  static  bending  response  of  statically 
indeterminate  structures.  A  boundary  sona  correction  approach  has  been 
presented  and  illustrated  which  permits  localised  boundary  restraint  to  be 
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accurately  modeled. 

kCMtuattd  none laaai cal  axial  atraaa  affacea  ara  observed  in  tha 
raaponaa  of  hyfrotharsally  conditioned  orthotropic  beams.  Other 
engineering  thaoriaa  fail  to  predict  thia  interesting  behavior. 

Baaed  upon  the  fiadiags  of  present  work)  there  ara  logical 
auggaatioaa  for  future  research. 

1.  The  study  of  alternative  clasping  definitions  has  indicated  tha 
presence  of  so  edge  sons  sad  decay  length.  These  are  auch  no re  pronounced 
in  orthotropic  structures.  A  reappraisal  of  composite  aster ials  testing 
aethods  when  the  speciaeas  are  relatively  short  is,  therefore,  a  worthwhile 
study. 

2.  Another  important  area  of  practical  interest  is  buckling  and 
poetbuckling  behavior  of  advanced  composite  structures.  A  large 
displacement  theory  based  upon  the  new  equations  is  recommended  to  be 
developed  to  pursue  the  above  issue. 

3.  Finite  element  models  based  on  the  present  theory  would  be 
useful  in  the  numerical  analysis  of  practical  structures. 
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APTUDIX  A 

OKTHOnoriC  SLAB  tOLOTIM 

! 

A  plana  at rear  aolution  for  flexural  «vi  propagation  in  an 
orthotropie  baas  ia  developed.  Thia,  by  auitabla  modification  of  alaatic 
cooatanta,  can  ba  uaad  to  obtain  atraaaaa  and  phase  velocity  in  an 
infinitely  vide  alab.  Tha  raaulta  ara  uaad  in  tha  dynamic  tbaory 
validation  atudy. 

Tha  looka'a  Lav  for  orthotropie  materials  v ay  ba  expressed  aa 


*11  J«,«  “•'» 

"»  *  "’xJ  *.,l  <*•» 

*  G1J  I",.  *  *,*  I  <*•» 

Tha  above  ara  obtained  by  invar  ting  Equation*  (27)  *  (29)  and  uaing 
Equetiona  (7).  and  ijj  ara  givan  by 

*11  *  *ll,(l  “  V13  V3l)  U,4) 

*33  “  *33,(l  "V31V13)  (A’5) 

Tha  aquationa  of  motion  in  two  diaanaiona  ara 
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Ut  tlM  b«M  vibrtt*  with  circular  frequency  w  .  The  loweet  eod« 
hee  one  ha  1  f  -wave  in  the  i  direction.  The  dieplaeenant  component  a  u  and  w 

nay  be  taken  ee 
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0|  and  »j  are  arbitrary  conetante.  p  ehould  be  chosen  each  that  the  above 
aatiefy  Iquetione  (116)  and  (117).  Thie  condition  leada  to 
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Ut  pj  Ml  pj  b«  tho  root*  of  IfitttiM  (A  .8) .  Tboo  tb«  aolutioo  for 
«  Ml  w  way  bo  upniMl  M 


«  •  cot  21  iia  ut  |*Aj  ciali  -j-  *Aj  tinb  j  (A.  13) 

¥  *  tin  2jr  ain  «tt  £  Aj  coch  -j-  ♦  A^  eoeh  j  14) 

A|  *  A^  art  Arbitrary  eooataata.  lha  rolatlooa  among  thoao  arc  obtained  aa 
follawat  tubatitotioa  of  I^MtlMa  (A .8)  and  (A. 7)  into  (117)  uaiog 
Bquatioaa  (A.  2)  am d  (A. 3)  load a  Co 
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equation  (A. 13)  ouat  bo  aatiafiod  iodopandantly  bp  Cho  oolutiooa  for  o  and 
«  aoaociatod  with  tho  root  a  Pj  and  p^.  Tbaraforo  fro*  (A. 13)  for  the  root 
pjt  it  ia  required 
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Bq nation  (A. 16)  ia  espraaaed  ia  Cho  followiat  term  toe  convenience  in  the 
aubaoqaaac  aoalpaia. 
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fr*qa*ocy  d*t*raia*at  by  m!«|  cb*  itrcii  booadary  cooditloo*  |i*w  la 
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